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ABSTRACT 


The extreipe sj ee 1 and altitude ranges of which modern^ high perfor¬ 
mance, supersonic air[ lanes are capable liave become manifested in wide 
variations in airframe transfer function parameters as flight conditions 
change* Airframe control systems and stability augmentation devices 
utilized In these aircraft require some form of gain changing to pro¬ 
vide satisfactory stability and control characteristisc over the complete 
flight envelope* In present day operational aircraft, system gains are 
scheduled as functions of airspeed and altitude, however these mechaniza¬ 
tions have seme shortcomings which are pointed out in this paper, and also, 
the desirability of self-adaptive aircraft control systems is discussed. 

A survey/ of adaptive techniques which have been applied to the airplane 
control problem is presented, and a self-adaptive pitch-stability aug- 
mentor which utilizes a digital computer In the adaptation loop is pro¬ 
posed, The proposed seIf-adaptive system is tested in conjunction with 
a digital three degree of freedom flight simulation of the pitch dynamics 
of a representative supersonic fighter-attack aircraft, and is found to be 
capable of performing the adaptation process in tv/o cycles of airplane 
short-period oscillation. Requirements for practical mechanization in 
a flight vehicle and recommendations for future research are discussed* 
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- UTTHOWCTIOU 

As the speed at: 1 c citude capabilities of modern high performance- 
supersonic aircraft have increased in recent years, the range of varia¬ 
tion of the dynamic characteristics of these aircraft have increased pro¬ 
portionately. A ten to one variation in damping ratio and natural fre¬ 
quency of some of the dominant modes of aircraft motion is not uncommon 
for modern airplanes with Mach 2 and 70,000 ft. altitude capabilities. 

As an example, for the airplane considered in this study, the damping 
ratio of the dominant longitudinal mode varies from ^ = 0.445 at M = 0.9 
0 10,000 ft. to^ = 0.064 at M = 1.2 f? 45,000 ft. 

Various stability augmentation devices have been used to decrease the 
range of variation of airplane parameters and to hold them as near as 
possible to optimum values. These devices, in their historical order of 
application, include fixed gain dampers, pilot-adjusted-gain dampers, 
speed-and-altitude-programmed-gain dampers, and seIf-adaptive-gain-adjust¬ 
ing dampers. A number of methods have been investigated for mechanizing 
self-adaptive stability augmentation systems for aircraft, utilizing vari¬ 
ous techniques which ’v^rill be discussed later. And, the purpose of this 
study is to Investigate the use of an airborne stored program digital 
computer for generation of the adaptation commands in a self-adaptive 
pitch damping system to be used in a high performance fighter-interceptor 
type aircraft. The dynamics of the airplane, including control system 
components, and the adaptation computer were simulated on a Control Data 
Corporation 1604 general purpose digital computer. 

Befoi^e beginning the discussion of adaptive control systems as applied 
to aircraft, it is felt that a better understanding of the system require¬ 
ments may be gained by a short review of the characteristic dynamic modes 
of a supersonic aircraft and the theory underlying the analysis and syn¬ 
thesis of these motions. 
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T: cr TV.£ 

Thf? ilx'plaue t. i-- axis system vhirh is generally used in airplane 
stabiHtv raid centre] studies is shovai in figure 1, as three inutually 
perpendicular axes intersecting at the airplane center of gravity. The 
longitudinal, or X, axis lies in the airplane's vertical plane of sym¬ 
metry and is atbitrarily oriented by the airplane designer so as to de¬ 
fine a horizontal plane of symmetry of the fuselage. In almost all air¬ 
craft the X body axis does not coincide \>;ith the longitudinal principal 
axis (longitudinal axis about which the moment of inertia is least) due 
to the high location of the mass of the tail surfaces and due to the 
assymetrlc top-to-hottom configuration of the fuselage. Consequently, the 
longitudinal principal axis normally slopes downward from tail to nose and 
coincides with the X body axis only at the center of gravity. The dynamic 
unbalance of mass about the X axis thus causes the airplane to have an 

1 product of inertia which will be important when cross coupling of air- 
XlU 

plane motions is discussed. 

The Z body axis, or normal axis, is perpendicular to the X axis and to 
the wing plane and has positive direction downwarvl. And, the Y, or lateral 
axis, is simply defined as being normal to the X and Z axes. Fortunately, 
due to the lateral symmetry, an airplane has no dynamic unbalance about the 
Y axis and thus only one product of inertia, exists. 

A second set of references, the earth reference system, must also be 
considered due to action of components of gravity along the three axes. 

To determine the gravity comp.cnents, the airplane's attitudes in pitch and 
roll must be determined and the acceleration of gravity resolved along the 
three axes by utilizing the Euler equations. The pitch and roll angles are 
illustrated In figure 2, v;here pitch angle is show with wings level and 
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AXIS 

X 

y 

Z 

Units 

1. Forces acting on the airplane 
along the axes: 

X 

x 

z 

lb. 

2. Moments ' acting on the airplane 
about the axes: 

L 

M 

N 

ft-lb. 

3., Moments of inertia of.the 
airplane about the axes: 

Ix 

ly 

'•Iz 

* \ 

slug-ft^ 

4. Linear velocities of the 
airplane's CG in space: 

u 

V 

w 

ft/sec 

5. Angular displacements of the 
airplane about the axes: 


e 


rad. 

6, Angular velocities of the 
airplane about the axes:.. 

P 

q . 

r 

rad/sec 


Figure 1. Body axis system and parameters used in the 
digital simulation 
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In tv;o ronflltlons the earth 


' -u-' .n 1.' (r^-y ’ev.-*) 

referetice anf,les *^.re. eqjlVt’Ient to Inte^r^tions of airplane angular vel¬ 
ocities about rheir r* spective boJy axes. But, if the craft is in an 
attitude with non-zero values of both pitch and roll, these angles must be 
obtained froni the projection of the X and Y axes onto two vertical planes 
containing these axes, 

A third earth reference angle which has not been illustrated and which 
plays no part in the dynamics of motion is the yaw, or heading, angle which 
is the projection of the longitudinal body axis onto the horizontal plane* 

Characteristic Airplane Dynamics 

The six degrees of freedom of ir.otion of an aircraft can now be defined 
as three angular freedoms in roll pitch and yaw about the X^Y, and Z axes 
respectively; and longitudinal, lateral, and normal translations along the 
X, Y, and Z axes respectively. Thus, the dynamics can be represented in 
terms of six simultaneous differential equations of motion. Before dis¬ 
cussing the makeup of the six equations of motion, it is felt that a better 
understanding of the following material may be obtained by a review of the 
various modes of aircraft motion. The motions can generally be divided in¬ 
to two classes, longitudinal modes and lateral-directional moder* The first 
class involves pitch rotation combined with translations along the X and Z 
axes, whereas the lateral-directional modes involve rolling, yawing, and 
translation along the Y axis. This division of modes occurs because of the 
low level of aerodynamic and Inertial cross coupling betv;een longitudinal 
and lateral motions. P'or small amplitude oscillations (excursions in pitch, 
roll, and yaw less than 15-20 degrees) the two classes can be assumed dynami¬ 
cally independent. 
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Figure 2. Earth reference system used In the digital simulation 
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* s 'C s }< ^-rr; |1or;s f’le \\ ,i ^ I' 1 

diracr lon.Tmo i^s ot o r:. rr^sei^tatlve ?=;uperr.onLc nlrplane,o The top haJi: 
of the ^ pl^ne Is shown for eacVi case. It is seen that the longitudinal 
motions occur in two oscillatory modeSj one with a relatively long period, 
the phugold mode, and one with a nnch shorter period, the short-period 
mode. The damping of the phugold is usually lov^ and in some cases it is 
a neutral or even divergent oscillation. However, due to its long period, 
the motion Is easily kept under pilot control with little or no conscious 
attempt at correction, and no artificial stability augmentation is re¬ 
quired for this mode. The phugoid occurs at coUvStant angle of attack vand 
manifests itself as a continuous trade-off between airspeed and alt tude. 

Figures 4. and 5. show that the period and damping ratio of the short 
period mode of a supersonic aircraft can vary over wide limits, and the in¬ 
clusion of the tlnie lag of the human pilot in the control loop can some¬ 
times cause the short-period poles to move into the right hand side of the 
S plane. This necessitates the use of pitch stability augmentation not 
ool)’' for improvement of handling qualities but for reasons of safety. 

This mode occurs at constant airspeed and appears as a pitching oscilla¬ 
tion accompanied by small amplitude translations up and down the Z axis. 

It is extremely important in maneuvering flight and is of prime interest to 
designers of airborne fire control and stability augmentation systems. 

The lateral-directional pole-r,cro plot for a supeisonic aircraft 
normally exhibits only one pair of ccmplex poles and tvc real poles as 
shown in fig. 3. The frequency and damping ratio of the lateral-direction¬ 
al short period, or Dutch roll, mode vary as flight conditions change in a 
manner similar to the variation of longitudinal short-period characteristics. 
This mode is a combined rolling-yawing cscillatic>n with small amj lltude 
lateral translations, ^nd its name Is derived from the swaying motions of 
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Figure 3 Pole-Zero Plot of the Characteristic Modes of 
Motion of a Supersonic Aircraft. 
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a Dutco skf'ter* The 3u!:cb roll oscilintions of a supersonic aircraft 
usually requure artificial darr^ping in the same manner as required by the 
longitudinal short-period mode. Therefore^ even though this study is con¬ 
cerned with an adaptive longitudinal stability augmentation system, the 
results and conclusions are applicable to the lateral-directional case 
as well. 

The roll damping, or roll subsidence, mode is shown as a negative 
real pole and it repi'esents the viscous roll damping of the aircraft. 

And, the spiral mode is usually a steady divergence which manifests itself 
as a gradually steepening turn with accompanying sideslip. The divergence 
is fortunately slew and, as in the case of the longitudinal phugoid mode, 
is easily controlled by the pilot. The roll subsidence and spiral mode 
poles can occasionally form a complex pole pair in the case of some sub¬ 
sonic airplanes v/ith large dihedral effect. This lateral-directional 
*^phugoid” mode will exhibit small amplitude, long period rolling oscillations 
with considerable lateral translation and very little yawing oscillation. 

This concludes the survey of normal modes of airplane motion. How¬ 
ever, there are other ^pathalogical” modes v^hich can be experienced when 
the rates and amplitudes of motion become large. Some of these are the 
phenomena of inertial cross coupling in roll, autorotation (Ref. 125.), 
and longitudinal pltch-up (Ref, 86); all of which can be eliminated or al¬ 
leviated by the use of well designed airplane stability augmentation systems. 
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Aerodynanic Stability Derivatives 
One final topic, the concept of aerodynamic stability derivatives^ 
should be mentioned before proceeding to the equations of motion. The 
forces and moments which act on an aircraft are generally expressed in 
terms of non-dimensional coefficients multiplied by appropriate aircraft 
and flight condition parameters. As an example, the equation for airplane 
lift is shoX'Tn below, 

P 2 

Lift = ^ VS 

iJhere: 

Lift is in pounds 

C is the airplane lift coefficient 

Li 

p =; atonispheric density in slugs/ff 

y St airplane forward velocity in ft/sec. 

2 

S airplane wing area in ft . 

P 2 

The expression ^ V is a measure of the kinetic energy of the air- 
stream and is termed the dynamic pressure because it also represents the 
force which would theoretically be exerted against a unit square surface 
held normal to the stream. Thus, the product of dynamic pressure and X7ing 
area represents a hypothetical aerodynamic force from which the actual air¬ 
plane lift can be obtained through multiplication by the lift coefficient. 
Moments acting about the aircraft can be obtained in the same manner with 
the inclusion of a moment arm, for dimension considerations, as is done in 
the following pitching moment equation. 

P 2 — 

Pitching Moment ^ C tt V Sc 
^ m2 

VJhere the moment arm c is the wing mean aerodynamic chord. Corres¬ 
ponding expressions for drag force, side force, rolling moment, and yaw¬ 
ing moment coefficients are included in the table of symbols. It is seen 
that there exists an aerodynamic coefficient for each of the six degrees 
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c>£ freedom of luotion. 


These coefficients and their derivatives are 


usually determined eiifpirical ly in wind tunnel and flight studies be¬ 
cause of the inaccuracy of analytical determinations. 

Each coefficient can be expressed in a Taylor^s series such as: 

. 

. 

The effects of the second order snd higher terms are insignificant in aero¬ 
dynamics and are dropped. So, for practical purposes the coefficient can 

, I ^Cl \ I I dCiS 


be expressed: {dQjJ\ ' ■ i vv^i-i / 


CL=te; a- \^l cc CJ • 


The partial derivative terms such as *1’^® called aerodynamic stability 

derivatives and are usually written for brevity in the form 
. . Ci= Ci^cc d- Ci_^ CO -f Ci^ 

Figure 6, a curve of lift coefficient vs. angle of attack, can be used to 
illustrate the physical meaning of stability derivatives, Here,(^ ^ - ) 
is the slope, or derivative, of the airplane lift curve* It is seen that 
the curve is linear up to CC 12"’. This indicates that - C CC is 

L hQr 

physically valid for small variations in (X • The fact that almost all 
aerodynamic stability derivatives remain nearly linear over a similar range 
of their corresponding variables enables the small amplitude motions of the 
aircraft to be expressed as a set of simultaneous linear differential equa¬ 
tions. 

How many stability derivatives must be considered in the equations of 

'k k 

motion'^ Given the six force and moment coefficients , C,,, ^ C|^, 


"^The lift coefficient and the drag coefficient C act perpendi¬ 
cular and parallel, respectively, to the relative wind, out it can be 
assumed with little error that they are. also normal and parallel to the 

X body axis. Then C_ = and 

Z D 
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0 . 


- 0 . 


- 0 , 


- 1 , 

-20 -10 0 10 20 • 

Angle of Attack, deg. 

Figure 6. Plot of Airplane Lift Coefficient vs. Angle of Attack, 
and and the fact that each coefficient can be"influenced by: 

1. angular displacement about each axis, 

2. angular rate about each axis, 

3. velocity along each axis, and *' 

4. acceleration along each axis 

there are 12 possible derivatives for each coefficient or 72 in all. Also, 
control surface deflections can add 36 more terms. -^Fortunately, airplane 
sjrmmetry, which causes a lack of cross coupling between longitudinal and 
lateral - directional modes, greatly reduces the number. The magnitudes of 
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some of the remaining deriv%atso small as to itake them insigni* 
fleant. As a result there are approximately 25 derivatives considered 
important enough to be included in a simulation program. These are surn-* 
marized in Appendix TIT. A list of stability derivatives for a representa¬ 
tive supersonic aircraft is tabulated in Appendix I. The aircraft is the 
same one for which the longitudinal short period mode characteristics are 
illustrated in figures 4, and 5. 

The Airplane aquations of Motion and the Longitudinal Short Period Transfer 
Function, 

The six equations of motion corresponding to the six degrees or free¬ 
dom of an aircraft are derived in appendix IV and the final results are 
presented in table I. It should be noted that the equations represent the 
balance of inertial and aerodynamic forces and moments along and about the 
airplane body axes, and that the terms are composed of airplane physical 
constants, aerodynamic stability derivatives, and aerodynamic variables. 
Assumptions used in the derivation of the equations are: 

1. The change in aerodynamic forces and moments due to deviations 
from an equilibrium condition are assumed directly proportion¬ 
al to the deviation. 

2. The deviations of the airplane from equilibrium are assumed 
small, 

3. The airplane is a rigid body. 

4. Constant thrust is assumed during an airplane distrubance, 
and secondary power effects are neglected. 

The equations are used directly in the digital computer simulation of the 
airplane dynamics, v;hich is outlined in appendix VI, and the airplane 
longitudinal short-period transfer function is derived from the pitching 
and normal translation equations of motion in appendix V. This transfer 






function is of prime impoi'tance in the study of the airplane pitch damper 
because the short-period mode is the dominant longitudinal mode* The 
short-period transfer function is written below both in its general form and 
also in terms of the airplane physical characteristics, airplane stability 
derivatives, and aerodynamic parameters. The second form of the transfer 
function Is shovTn to Illustrate the dependence of the dynamics of the short 
period mode on changes in stability derivatives and aerodynamic parameters, 
and it is possible from inspection of the function to determine which de¬ 
rivatives and parameters exert influences on a given characteristic, such 
as damping ratio, of the short-period mode. 


Longitudinal Short Period Transfer 


Function 


A. General Form 




e _ Ke(5^Te) 


B. Detailed Form 
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sin^ cos & , ^ pitch angle 

50 roll angle 

coBtp cos Q , j 1 

' (Angles between the body axis 

and earth gravity axle) 
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The ptsr|>ose of ^^Itplane st^tblllty augnuentiti ion is to reluce the 
range of variation of a iynafnic characteristic, or characteristics, and 
to hold the variables of interest ns close as possible to predeter^uined 
optirpurn values* In the case of an aircraft these optimum values are de¬ 
fined in terms of ’^airplane handling qualities’* which describe an air¬ 
plane’s response to manual control by a pilot. Even though the determina¬ 
tion of desirable handling qualities is subjective, and is dependent on 
pilot opinion, a number of studies have been conducted to specify these 
qualities in quantitative engineering tennis* One of the better known 
studies, conducted at the Cornell Aeronautical laboratory, (Ref* 34), con¬ 
cludes that the desirable natural frequency and damping ratio of the short- 
period mode of a fighter airplane lie in the range ~ to .5 cycles per 
second and damping ratio = .65 to .70. The same study also concludes that 
even though the longitudinal natural frequency is of importance In maneuver¬ 
ing flight, it can vary over wide limits and the handling qualities V7ill 
still be considered tolerable, as long as the short period damping ratio is 
held in the .65 to .70 range* Thus It Is considered more desirable to equip 
a high performance airplane with artificial damping than to provide static 
stability augmentation. Tt should be noted that for the craft considered 
in this study, as well as for almost all other airplanes of similar type 
and mission, the longitudinal damping ratio of the un-augmented airplane 
is always belov/ the desired .7, as is sho\^ in figure 4, This means 
that the pitch damping system will work to supplement the basic aerodynamic 
clamping in all cases. 
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Mechanization oC the I itch harper 


In the discussion of the airplane longituUnal sliort-period transfer 
function it vas sho\^ that the short-i)eriod motion can be represented as a 
second order system v;ith characteristic damping ratio and natural frequency. 
It was also shown that the natural frequency and damping ratio can be ex« 
pressed «is functions of the aerodynamic stability derivatives. The ex¬ 
pression for short period damping ratio being: 


A r - r 
ms/ ^Loc. avrv^^q^cr 



From the above relationship it can be seen that the most effective means 
of varying the damping ratio is to artificially vary the pitch damping de¬ 



resistance to pitching velocity. 

The pitch damper must then provide an aerodynamic pitching moment which 
is proportional to airpLane pitch rate, and v^bich has a direction in opposi¬ 
tion to the angular rate. Thus, the moment produced is in phase with the 
viscous aerodynamic damping in pitch and the tv;o moments combine to form 
the total pitch damping moment. The artificial damping moioent Is usually 
generated by deflecting the airplane® s horizontal st^’tbilizer in direct 
proportion to measured pitch rate. 

Figure 7. is a block diagram of the pitch dampe'^ u^t.-rd In this study^ 

The airplane transfer f'metion shov;n In the diagram represents the longi¬ 
tudinal short-period mode only because, as will be demonstrated later, the 
phugoid mode and the lateral-directional dynaralcs are of negligible impor¬ 
tance in longitudinal stability augmentation considerations. 

The stabilizer actuator in the diagram is an iri^averslble hydraulic position 
servo, and the parameters in the t ransfer function rev resent the state-of- 
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Figure 7« Block Dlt^.grarri of the Pitch 
the-art €is used in modern aircraft. It is generally assumed that the 
dynamics of the rate gyre have negligible effect on system operation be¬ 
cause the ratio of the characteristic frequency of: stale-of-the-ait rate 
gyros to the short-period mode frequency is of the order of 20 to 30 
(Ref. 43), and tlms the gyro transfer function is assumed to be unify. 

Table II lists the 11 airplane transfer Functions use! in the study a Toni: 
with corresponding flight conditions. And, figures 8 thru 18 are root 
locus plots showing tlie behavior of the syste m vjith vi^rylng :! in at 

each condition. The square symbols on the plots indicate the galr values 
at which the airplane short-period damping ratio has been increased to 0.7. 
It is seen that the use of the damper has V3cy little influence on the 
short period frequency when operate! over the range of gains require.I to 
bring ^ up to 0.7 at each flight condition. And, it is seen that 


























the 


i^ole inoves ir 3i=ch a direction, with iwcteaching datoper gain^ 
that the actuator d^AViiivg ratio rer'alas close to 0*7 over the danper opera¬ 
ting range. The actuator t»ode will become unstable at extremely high gains, 
but operation at such gains^ is of no interest In this study. 

It is knov;n that the closing of a feedback loop around a plant having 
varying characteristics will cause the plant output to vary over a smaller 
range, as plant characteristics vary, than the range of variation of out¬ 
put ivith no feedback. From this it x^ould be expected that the pitch damper 
should not only Increase the average value of damping ratio over all flight 
conditions, but also that the range of variation of damping ratio, as condi¬ 
tions vary, should decrease. This is borne out in Figure 19. x/hich is a 
plot of short-period damping ratio vs. speed and altitude with constant 
damper gain s 0.3. The chosen value of damper gain is approximately the 
median of the gains required for 0.7 damping ratio over all 11 flight 
conditions. In the case of the basic airplane, ^ varies frorti 0.445 at M 
= 0.9 (? 10,000 ft. to ^ = 0.064 at M s 1.2'^ 45,000 ft. whereas the varia¬ 
tion of / uith damper gain = 0.3 Is from l.l to 0.436 for the same two 
flight conditions. Thus the range of variation about the mean has been 
reduced from approximately 190% to approximately 95%. 

Even though the use of the constant gain damper has caused a consider¬ 
able improvement in the dynamics of the airplane, a further restriction in 
the range of variation of damping ratio is usually desired for fighter-inter- 
ceptor type aircraft. One method of solving this problem^ x^hich is in 
general use today, is the orogramining of stability augmentation system 
gains, as functions of measured flight parameters, in such a manner that 
the programmed gains will cause near optimum system response over all 
flight conditions. The techniques of gain scheduling methods range from 
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variaticn ol' JSfnper gain ar function of f.aasuned airspeed to systems 
which program as functions of airspeed* altitude^ fuel load^ center of 
gravity position, wing flap position, landing gear position, and vari^ 
cus combinations thereof* One fact becomes obvious from the preceeding 
discussion; the dynamic characteristics of the aircraft must be accur-- 
ateiy knov;n in advance at all flight conditions and for all airplane 
confLgui*ations if the progranuning system is to operate successfully 
Such a detailed determination of airplane characteristics requires a 
large amount of wind tunnel and flight testing vjhich is both time consuming 
and costly* And, even after the desired information has been determined 
it is usually necessary to conduct lengthy flight tests v;ich the program¬ 
ming and stability augmentation systems installed in the aiirplane for final 
’’tailoring'* of »::he aircraft-system combination* In addition, accurate in¬ 
flight air data and airplane configuration Infoni^ation must be continusouly 
fed into the system during flight to provide the basis for gain scheduling* 
The instrumentation and transducers required are generally costly, and are 
usually complex enough that they create reliability and fail-safety problems. 

A possible solution to these control problems is to design a system 
which is intelligent enough to detemine when the airplanedynamic per¬ 
formance has changed, and In what manner it has changed, and to use this 
information to adjust the stability augmentation system gains to obtain 
optimum performance over all flight conditions. This is the concept of 
the self-adaptive control system. 
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IV 


SELF-ADAPTIVE CONTROL SYSTEMS 


The concept of adaptive control has been applied to a wide range of 
control system types with the result that a number of definitions of adap¬ 
tive control exist. Hov/ever, the most often used definition^ and also the 
one which appears to best fit the aircraft control problem is as follows: 

An adaptive system is a control system configuration in which 
the measurement of process dynamics or signal characteristics 
Is utilized to automatically adjust the controller in an attempt 
to achieve optimum operation at all times. 

From this definition it can be seen that the adaptive process must per¬ 
form two functions, the first of which is the ‘identification*’ function, 
or the measurement of process dynamics or signal characteristics to deter¬ 
mine present system state and compare it with the desired state. And, the 
second is the ’’command” function which is automatic changing of the system 
sc that it x^?ill approach a desired state. In most adaptive control studies 
it is assumed that the variation of system characteristics takes place slow 
ly vjith respect to the frequency of the system mode of interest, and thus 
the identification process can be allovjed a relatively long measurement in¬ 
terval during which to determine system state. Hov7ever, in the aii craft 
problem, it is desirable to perform the adaptation process relatively rapid 
ly In comparison to the quantity being controlled, in this case the short 
period mode dynamics. This is because virtually instantaneous changes can 
take place in the airframe transfer function due to release of payload or 
certain Mach number effects on stability derivatives. For this reason, 
one of the objectives of this study is the investigation of an adaptive 
technique capable of rapid identification and command. Another considera¬ 
tion in mechanization of the adaptive process is the generally accepted 
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reqaiiei«:ent t’u- ^ o :lei :oi sycteiM f/c he conSi.<levt: 1 ttaly adrip*^ 

tJ.ve, t;»e I'!entIficaiion yi*ocess c^:?j net involve tneasurencnl of any 
variab'^e which Is external to the sy^teiv. To Illustrate the co)}ce|,ts 
and requirei-ent3 discuvssei above, a block diagiam of Gibsonideal 
generalised adaptive system (Ref. 76) is shovm in Figure 2Cu 


IJentification and Como:.nd 
Computer 



Figure 20, Gibson Vs Ideal Generalized vdaptive 

System 

Classification of Adaptive Systems 
Tn ot'clei' to compare the approaches which have been taken to the 
design of adaptive control systems, Aseltlne, in reference has seperat 
ed them Into 5 classes which, in some cases, have a certain amount of 
overlapping, but even so, they form a convenient framework for comparison 
Aseltina's classes are: 
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1) CzSP^ive J:af:ior4: Fysteris uhicr. uitbout: autorp^ttc 

changes ii» sysreri! cluaracceristIcs, hut rather through de¬ 
sign fox' vinffono oreratlon over x^/lde variations In environ-* 
TTient. 

2) Input Signal Adaptation: Systems V7hich automatically ad¬ 
just their parameters in accordance with input signal 
characterlst ics 

3) Extremum Adaptation: Systems which adjust to obtain a inaxi- 
mum or ninirnum o£ sonje system variable, 

4) .System ^7artable Adaptation: 55y sterns irt xvhlch the command 
pi'ocess is actuated directly by measuremenf: of a system 
variable or variables* 

5) System Characteristic Adaptation: Systems In which the 
Identificat:ion process measures S 3 rstem variables in order 
to determine system! characteristics. The system character¬ 
istics thus deteimiined are compared with desired character¬ 
istics and tlie differences used to generate coranand signals. 

It has been noted from a survey of adaptive systems which have been investi¬ 
gated for possible use in aircraft control that^ to date^ all such systems 
fit into classes 4, or 5. And, it has also been noted that an alter¬ 
nate classification system might be applied to aircraft adexptive control 
systems. The suggested classes are: 

1) High Gain Adaptation: Systems in x^hlch the characteristIcs 
of the aircraft are submerged by the use of high feedback 
and/or forv;arrl loop gains. Some fornr of Idealized model is 
normally used to provide either system input signals or out¬ 
put comparison signals, and the techniques usually involve 
either nonlinear processes, with attendant limit cyclcj or 



i 





M 





operation of a n»lnor system mode near instabilityc 

2) Impulse Response Adaptation; Systems in which the lQenti'=* 
flcation process computes a time history of the plant Im¬ 
pulse response. The Impulse response Is then utilized by 
the command process to generate plant alteration signals. 

3) Faraineter Adaptation; Systems In which plant parameters 

are determined by direct computationo The error between 

the desired value of the parameter and its actual value is 

used as the command signal for plant alteration® 

Survey of Adaptive Techniques which hv"^ve Been Applied to the Airplane 

Stability Augmentation Problem 

1) High Gain Adaptation Systems 

The earliest studies of aircraft adaptive control involved systems 
which operated at high gains In order to force the variation of airframe 
parameters into a small rai^ge. One of the first of these studies was 

conducted by Rath of WADC (Ref. 134) on a method x^^hich was not truly ad¬ 

aptive, in that parameters were not automatically varied^ however, it 
may be classed as p^assive adaptation. Rath® s, system used a parallel com¬ 
plementary optimum response model as shox^n in Figure 21® The goal of 
this study was to design a system which would produce optimum airframe 
dynamic characteristics, as airframe parameters varied, with constant 
system gain. A fourth order transfer function was selected as the opti¬ 
mum model for the pitch dynamics and the same input signal was applied to 
the airframe and to the model. The difference between model and plant 
output was obtained and various functions of the difference signal were 
utilized to improve the airframe response. Both the forward loop gain 
and the shaping loop gain were made high in order to ’’force’* the system 
































Figure 21. Rath^s System for passive Adaptation Using a Corf;p lament ary 

Optimum Respons*^ Model 

overall characteristics to resemble the model. Feedback functions studied 
included pitch angle^ or position^ error Ceedbackj Integr^^ted position 
error feedback^ rate error feedback, and rate-error'-plus-positlorh-error 
feedback. It should be noted that error, as referred to here, means the 
difference between plant output and model output rather than pl^nt input- 
output error. The first two feedback functions produced overall instabik 
ity at high gains, whereas rate ertor feedback resulted In ^ stable system 
with Infinite steady state position error. Rate-error"plus-posiilon*arror 
feedback allowed very high gains to be used before the system became un¬ 
stable, and v;hile the system was stable, the airframe output rratcbed the 
model output acrept.^bly at most flight conditions while using a cotAStant 
value of feedback gain. It was also found that, whatev^ar feedback shaping 
was used, the basic airframe and hydraulic i?ctuator must fom a stable 
system with the shaped feedback loop alc-ne. That is, the f^3;e•'^f^L^rward loop 
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civjl*] rot. I'( of rhe ‘rfi -s.-. ’ sb.ap-,.i fV^'l^acl-: 

comhinatIcn* The optimufn i^fer furif.tlou feedforv^iioi it* comMnatlor 
v;itb tlui shape.! feeibcsc!* h*'d the effect cf carng^'n^ ihput to the 

alrfrai'e to a which c-?used Lbc airfran'e to have the desir:t.rl out¬ 

put, assuming tbs airframe hoc th^ desire! corstrcl rates and high gain 
control pcver« 

An approach sli-llar co "‘ath^s passive adaptation .method was taken by 
Jandols of Convair who utili^jed the inverse of the optimum response model 
in the feedback loop .In combination with a high forv/ard loop gain system 
as shown in Figure 22e (Ref* 43). 


AIRFRAME SHORT PERIOD 
TRANSFER FUNCTION 



Figure 22* Dandois’ oystetr. Using the Inverse of the Optimum >euSponse Model 

in the Feedback Loop 

He demonStrated that, with such a feedback arrangement^ the overall system 
characteristics approached the characteristics of the opciiium response 
model as forward loop gain was increased* In oper'ition^ the 

system held airframe para*>iters wichi^ acceptable Hr. its c -/iHe 
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r.5nge oi JHght coniiticns, nxid i .ghe-i ci.e .^nln »be Hetf er the 

perfornance* however, gain Krnlts, dependent upon "light conditions, 
were reached at which tlie systen becat^e unstable* The Inscability was 
shown to be due to sensor nonlinearities rather than a basic, theoretical, 
systefo instability. 

One of the first **high gaiiV^ adaptive systeiris to actually use an 
identification and cormnand process was the system-variable^adaptation 
system of Fluggc-Lotr. and Taylor (Fef. 64), which is illustrated in Fig¬ 
ure 23. As is shown in the block diagram, one of a number of alternate 



Figure 23. Nonlinear Adaptive System of Flugge-Lot:^ and Taylor 
feedback paths was used according to a switching logic process based on 
simultaneous measurements of system output, error sigtmlj and derivatives 
of these. The IdentlfIcaticn process thus consisted of the ^oeasurepants 
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S'*-' tch’U. •' '' il ••t'’ r>'• .-ipj. tb snd 

the COTcon*v-Gtte of Tv^.-. itching to the chosen path, Tcx^iird gain was 
kept high in order to ’*force” a small instantaneous error between input 
and output^ and the system v;as nonlinear during any length of time com-* 
parable to the natural period of the plant because of the switching be* 
tween feedbacl; loops during one cycle of operatlonp The arrangement was 
shown to be insensftlve to variations in plant paravfeters over wide limits, 
how*ever a limit cycle persisted during operation due to the feedback non¬ 
linearity. 

The Minneapolls*Honeyv7ell relay actuated adaptive system (Ref* 165)j, 
which is presently undergoing flight tests, is similar to the Flugge^hot 2 
mechanization but uses only one feedback path vith switching logic in the 
forward loop. The sw^itching logic derives Its input from the error signal 
between optimum response model output and the process feedback signal as 
shown in Figure 24. The switching logic actuates a relay^ or on-off con* 
troller, to operate the process at a high gain with the objective of 
maintaining a small error signal between model output and plant output. 

As in the Flugge-Lots system, the airframe output approximately dupll* 
cates the model output when the forward gain Is high^ and the character* 
istic limit cycle is also present. A second adaptive feature of this 
system, aside from the switching logic, is the use of limit cycle ampli* 
tude measurement tc vary forward loop gain. That is, forward gain is ad¬ 
justed in such a manner that the limit cycle amplitude is held constant 
over all flight conditions. The prototype of this system was success fully 
tested in an F-94C fighter aircraft and an improved version is presently 
being flight tested in the X-15 high altitude research airplane. 

Another application of th^ utilization of the inverse optimum re¬ 
sponse transfer function ixi the feedback loop, in conjxinction with 



Limit Cyo/e 


Amphtude. Mecx.sure>in^^^t 



ACTUATOR 


1 jRATE 


Igy/^o 



Figure 24. Minne^polis-Moueywe! 1 Fel^y Actu^'ted /^^.ptlA^e Systen^ 
high forward loop gain^, is the General Electric system (Ref. 32), vhich 
was recently tested in an F‘-102 interceptor aircraft. The mechanization 
is similar to that of Dandois shovii in Figure 22, except that an active 
system-variable-^adaptation fe/^^lure has been added. This f^^ature lnv*^lves 
variation of the forward loop gain, as flight conditions change^ to 
maintain the system just below the boundary of Instability. As was prt,^ 
viously meriticned, when consid^.rlng systems of this type^ the higher the 
gain, the closer i7lll be the airframe charactetistic.s to the optimuu^ re¬ 
sponse characceristics and thus if the system Is opernted ne^r the in¬ 
stability boundar}’, the tiaximum p*^rfon'^anie will he extracted. In oporA-' 
tion, the adaj: t&tion stem T.jonitors the frequency the control surface 
actuator because the actuator mode is the one which becomes unsta'* le as 
forward gain is increased. Tlse is then arlj isted to the 
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freqrt^incy of this r oc'e ccustr.iif over Flight cofMUtlons, It shoubi 

be riotea that the actuator mode Is a lightly dauiped osclllatloi and 
that non.ifil coritro! Inputs are utilf^.ed to provide the 'required pertur¬ 
bations^ The frequency rr.Gas<aemenl Is accomplished with a narro^v band 
pass filter and ^ero crossing counter arrangement« 

2) Impulse Fesponse Adaptation Systems 

The function of the identification process In the impulse response 
adaptation class of systems is first, the determination of a time Vtistory 
of the plant impulse response through measurement of system input and out¬ 
put signals over a period of timej and second, the detenmination of plant 
characteristics through examination of the Impulse response* An e«arly 
method of impulse response self optimization which was investigated by 
Aseltine (Ref, 4) serves to illustrate the process, ev^en though this 
particular mechanization would not be suitable for aeronautical use. 

The impulse response was obtalr.ed by periodic input pulsCwS which v/ere 
widely spaced with respect to the period of the plant mode of interesSt. 
After each perturbation, the system damping ratio was ccinputed by separ¬ 
ate integration of the positive and negative portions of the inpulse re¬ 
sponse and obtaining the ratio of positive to negative area. Damping 
ratio was cenipuLed from this ratio andi the system damping coefficient v;8S 
varied in a closed loop manner to cause actual damping ratio to approach 
a <leslred value. This ipethc 1 was considered undesirable for aircraft use 
prirr;arlly because nom-aT operating signals must be excluded frotr» the In¬ 
put during the comparatively lengthy identifIcatior period, 

A method of impulse response determination in which normal Input 
signals could be Incluiled In the tdentifIcntlon process w^vS Investigated 



bv /noerj^cn 
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i.?ci?^lr;u oF the 


convolution ’^beoxeur '.-beie’^v the r?y5^tev Inprtt convolved In tin}': v;ith tVie 
Impulstj tyspoure ixi ?, tme history of the system outfut^ In o^ete- 

t'ion, rati’oi!i noise with cjuiot carel cva tuinrt: on was suporlmyOv^ed ot\ 

the nocTr^al Input *15 shown In Figure 25* The plant output was £ilt- 

erej l:o sepArace the noirve co‘<'jnoi"jent fron» the tola* vesi and the f-ti- 

pat and output noise corpvuents were cross-correlated to obt’^ln ^)oil')- 

CO 

• The nets-? tiipit had t bandv?rdth 3 to 10 times the p? Tit 
bandwidth so that (pn(T) wis effectively an impu'*se hUa ^pioCT') Ar- 


proximtely the plant Iriipulse response* Twelve charniels of eorss ro>/rela-' 
lion \veve used and cl> cm^nnel provided one point on e time history of 
the imptilse response, Da.'ipi'ag >:atio was deterrdnel an I utilized as Jis- 
cussed in the preceeding pair^.grapli. 


PARAMETER PARAMETER 



Figure 25. 


Block riagra^'-i of ^ercot^tronic*s Impulse F.esporse ^eli Cvtlir.i. ' 
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':noLber rnelho - oc se jet*, rmui^icjon i ti.liic.ing the convolu¬ 

tion thaorer.T In conjunction wltli si/».ultaiieous iVifut ntnl output mennurenentf> 
Is currently under investigation at Hughes Aircraft e In this rriethcdj dis* 
Crete the input and output tlnie histories are use! as coeffici¬ 

ents in two polynomials. The output polynomial is periodically divided by 
the input polynomial thi'ough the use of a high speed digital coijputer. The 
coefficients of the resultant polynomial can be shouni to be amplitude sam¬ 
ples of the system impulse response. One requlierrient for this method is 
that the amplitudes of the systen output samples nust be corrected foz“ 
initial conditions existing at the beginning of the Interval over which 
samples are taken. The correction process requires a knowledge of the 
impulse response and thus the impulse response must be determined in a con¬ 
vergent sticcesslon of polynoml.nl divlstovis* Air present, the requirements 
for convergence to the true inpuise response are under investigation, 

3) Parameter Adaptation Systems 

The function of the identification process in a parametar adaptation 
system is the direct determination of plant parameters through measurements 
of, and operations upon, selected signals and their derivatives. 

One of the first studies of parameter adiaptaticn vjas conducted by Kal¬ 
man, of Columbia University, who used a digital computer to determine plant 
parameters (Ref, 91), In Kalman*s system, coefficients of the plant*s pulse 
transfer function were computed by continuously measuring and storing 
sampled values of plant input and output. After a series of samples had 
been stored, values of the transfer function coefficients were assumed and 
a series of sample output points were digitally calculated using the known 
past input samples. The mean square error between calculated and known out¬ 
put points vjas determined and the assumed values of the coefficients were 

'»6 



m 



alter-' * rl o 1s-ast souu r.s f il s.V:_. which 

cl?,angea ce^eificlents lii r* pv.aiiiner that terhKed the nean error'. This 

process was repeat'.-' At ver}.- high s| eed until the '‘rtE was I'educed belw a 
desired limit. After the f-lant coefficients had been determined with suf¬ 
ficient accuracy, an optimal controllei: was synthesized In a second digital 
computer routine arnl the results were used to adjust the system controller* 
A different approach to parameter adaptation was taken by Whitaker, of 
MIT, who used a model reference as shown iu Figure 26, (Ref, 171)• 


MEA6URED OUTPUT QUANTITIES 
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Figure 26. MIX Model Reference Parameter A-Japtatlon System 
The input signal was fed to both the plant and f,o the Ideal model, and a 
difference function between plant output and model output was used to form 
an integral squared error. The ISE thus formed v'as required •■o be related 
to a single jilant parameter, which was then varied in a closed loop manner 
to cause the tSE to attain its niinimum value. The primary objective of the 
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study was the deternri.nrjtilon of vjhat functions of the plant-output-to-model~ 


output error were Indicative of what parameter variations. After the de¬ 
sired relationships were determined, the system was mechnnir.ecl and fll;;;ht 
tested in an F-94.\ fighter aircraft and it was found that adnatation could 
occur within an interval corresponding to two or three time constants of the 
dominant system response modes. However, recent developments have made 
possible a significant reduction In convergence time and the elimination 
of much of the complexity of the previous mechanization (Ref. 127). These 
developments were derived fircm a recent doctoral research y rograrn at -^IT 
which evolved an analytical design method for parameter adaptation model 
I'eference systems. 

k iriethod o£ adaptation etuploylng an analog computer, which calculated 
values of tu-tantaneous for^^rrd loop gain, was stiidled by Corhln of Broo! - 
lyn Tolytechnic Institute (Rtf. 38)* The gain value calculated was com¬ 
pared with n de^sired level and the difference used as a signal to the gain 
controller* The differential equacioris of motion of the plant were set up 
in the analog comruter as Integral equations and were solved for the gain 
value, given continuous measurements of the plant input and output signals, 
as shex^ In Figure 27, 



Figure 27. Corbin* r Gain Ada|tation Method TTtili^ing an Analog Computer 

in the Identlfleal Ion Irocess 























A Miethod simlUix' lo Cox’bU/s npproach, but x^hich eniployed a digital 
computer to directly solve for the plant parameters, v;as investigated by 
Azgapetian of Servomeclxanisrns Incorporated (Ref, 7*)* 1^ this study, the 

equations of motion of the second order system under investigation v^ere 
converted to sample data notation, and relationships for damping ratio, 
natural frequency, and forward loop gain were derived in terms of past 
sample points of plant input and output« The equations i^ere then mechaniz¬ 
ed in a digital computer routine which periodically sampled the input and 
output signals and coiuputed plant parameters at a high rate* parameters 
thus calculated vjere used in a closed loop manner to vary a controller which 
caused overall systetA dynamics to remain constant as plant characteristics 
were varied. 

The adaptation method investigated in this study is similar to that 
discussed in the preceeding paragraph in that a digital computer is util¬ 
ized to solve second order equations of notion for system damping ratio 
and natural frequency given system output samples. The primary differences 
are the method of sobition, and the fact that the system under investigation 
is actually a sixth order system which includes a second order hydraulic 
servo and a fourth order airframe transfer function. 
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V - i^OTO^^D SYSTEM 

As mentioned In Che preceding section, It is proposed that an air¬ 
craft pitch damping system could be made to produce invariant short period 
damping ratio over all flight conditions by digital computation of the 
short-period d«smplng ratio and closed looj) adjustment of the pitch damper 
gain to cause the damping ratio to equal a desired value* A block diagram 
of the proposed system is shown In Figure 28 below* The proposed identifi¬ 
cation process thus consists of solving the airframe equations of motion 



Figure 28. Freposed Self-Adaptive Fitch Damper 
Cor values of damping ratio given samples of pitch angle and its deriva¬ 
tives, and the command process consists of varying datnper gain at a rate 
proportional to Che error between desired and computed / and in the 
direction indicated by the sign of the error* Basic assumptions and condi¬ 
tions for this study are: 
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1, The Alrfrpiiie and control system are linear. This has 
been sbovn to be true for small perturbations* 

2. There Is no control surface rate limiting or deflection 
limiting. 

3. Airframe parameters are time variant but are assumed con¬ 
stant over the length of time required for one identifica¬ 
tion cycle, 

4, System disturbances due to external noise (atmospheric 
turbulence) will not be considered* 

It is also assumed that the short-period mode is dominant in the longitudin¬ 
al dynamics and thus the airframe longitudinal response can be closely des¬ 
cribed b}" the short-period frequency and damping ratio. This assumption 
V7as verified by conducting two longitudinal flight simulation programs at 
each of the 11 conditions listed in Table 11, The first program was run 
In two degrees of freedom, pitch and normal translation, and thus simulated 
the short-period mode dynamics, and tlie second program simulated the short- 
period and phugoid modes by inclusion of the longitudinal translation degree 
of freedom, Dansping ratio and period v/ere deternlned from pitch angle time 
histories obtained by pulsing the aircraft about the pitch axis and alloving 
it to return to ecuillbrium. Characteristics obtained from the simulation 
programs are listed in Table III, It is seen that the two degree and three 
degree of freedom characteristics are in close agreement, and thus the deci¬ 
sion was made to investigate a computational method of determining the 
■■pseudo** second order damping ratio of a sixth order, three degree of free¬ 
dom, airframe-control system combination from sampled data. 
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Table III - Sbovt-l ex j.oJ lio’e Chaj;actGt;Lsclea Comi arel with Dntr.j ing K?.tto 
arc’ Icriocl Obtained froi.. a bongitiidina 1 Three Degree of Free- 
dora Flight Sinnlation, 


Flight 

CoTid'!’t LOii 

Jti , 1 ^. irif. 

Ra*: to 

Pei^iod, 


M 

A1 1 

2 DF 

3DF 

2DF 

3DF 

0,3 

I0,CCC^ 

0.444 

0.440 

1.650 

1.700 

1.2 

10,000' 

0.185 

C.177 

1.050 

’ . 005 

1.4 

lo.ocn’'- 

C. 155 

0.156 

1.025 

1.020 

0,9 

25,non' 

0.245 

0.244 

1.974 

2.025 

1,2 

25,000" 

0.100 

C.1028 

1.225 

1.275 

1.4 

25,000" 

0.104 

1.105 

1.199 

1.175 

l.S 

25,000" 

0.095 

0.096 

1.200 

1.225 

0.9 

45,000" 

0.140 

0.141 

2.900 

2.075 

^ ? 
t. . 

45,000^ 

W , i 0*4 

0.066 

1.725 

1.725 

1.4 

45,000^ 

0.088 

0.067 

1.650 

1.675 

1.8 

45,000^ 

0.057 

0.C58 

1,600 

1.599 
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The generalized differential equation of motion for a second order 
system is: 

e^^^o)nB^w,^e= 

from which an expression for damping ratio can be derived as; 

/= Aife _ Q _ cor,e 
^ 2.u>ne ^uj„e ze 

It is seen that, in this form, the damping ratio is a function of the in¬ 
stantaneous values of 6 and its first and second derivatives, the input 
forcing function, and the natural frequency. Azgapetian has shown in Ref. 
7. that an expression for, / is independent of (A)n , can be de¬ 

rived as an integral equation. Bat it was decided, for this study, to 
avoid long term integrations and tr;eii attendant drifts, and rather to make 
use of the high speed iterating capabilities of digital computers by de¬ 
riving a companion expression for natural frequency to he used in an itera¬ 
tive process for arriving at ^ and simultaneously. The expression 

for natural frequency being; 


S ~ 


W 4 ‘- + 




KSe 

e 


which is In terms of instantaneous values of 0 and its first two deriva¬ 
tives, control input, and damping ratio. It is noted that the term, k, 
by which the control deflection is multiplied, corresponds to the airplane 
elevator effectiveness, or steady state airframe gain, which is also a func¬ 
tion of flight conditions and which therefore must be determined in a 
spearate computation routine. Corbin, In Ref. 38. has devised a method 
of computing the instantaneous value of forward loop gain of a second order 
system by utilizing integral equations. Corbin”s system utilized an ana¬ 
log computer which was supplied continuous values of plant input and out¬ 
put, but, the method could be converted to digital form with sampled in¬ 
puts and used in this study if desired. It was decided not to introduce 
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this additional coTmleKli"' Irto the rro;:;ran l^nt i-nthcr to concentrate 
solely on t'le plant parameter computation method and its application 
to an adapt ii'e system because the coraputation method was, as yet, untried. 
This decision necessitated the enclusicn of pilot control inputs to the 
system because, if pilot inputs were to be included, the value of the air¬ 
frame gain would have to be known. It should be made clear that control 
deflections called for by the pitch damper are not included in the preceedin 
expressions for natural frequency and damping ratio because the airplane 
pitching moment created by the damping system is considered to be added to 
the natural aerodynamic damping moment to make up total damping. If pilot 
inputs were to be included, the control deflection created by the pitch 
damper would be subtracted from total deflection and the remaining control 
deflection would be used as the input term, ^ , in the / and equa¬ 

tions, With pilot control input excluded, the equations are reduced to: 

f 

_ __e _ uiqQ 

ZUnG 2.e 

The proposed iterative computational routine is as follows: 

1) Measure and store sampled values of 0, G, and 0 at tv?o 
separate instants in time, 

2) Assume a value of damping ratio and compute ^r\ using the 
first sample set of 0, 0, & 0 and assumed . The 61^7 
equation is used in the initial calculation because the 
value of COf) is dependent upon variations in ^ to a lesser 
extefit than is the value of ^ dependent upon variations in 

(Mr, o 

3) Use the value of computed in 2) and the second sample set 
• • # 

of 0, 6, & 0 to compute ^ . 
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4) 'Tse the no:'' value of ^ and tie firat sample set to com¬ 
mute a new value of 00 . 

5) Repent steps 3) and 4) in sequence until the successive values 
of damping ratio and natural frequency either converge or 
diverge. 

Although it might seem that only one sample set of 0 and its first tv;o 
derivatives is required to define the parameters of a second order system, 
it is seen that the method described above requires two sets for convergence. 
If only one set v^ere used in both equations, the obtained in the first 
calculation, when substituted into the second equation would produce the 
originally assumed ^ , and convergence would be impossible. This is 
because the relationships for ^<7 and ^ are actually just two forms of 
the same equation. 

An appropriate question to ask at this time is, will the computations 
converge, and if so, under what conditions will they converge? The answer 
is given in the following section. 
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VI APPLICATION’ OF PPOPOSED TTSPATIVE DIGITAL DAMl^TNG PATIO COMPUTER TO 

A SECOND ORDER SYSTEM 

It was felt that it would he better to test the proposed plant identi¬ 
fication method first witVt a simple 2 nd order system before adding the com¬ 
plexity of the phugoid and actuator modes of the three degree of freedom 
flight simul 4 ^tion. When done in this manneTj the validity of the basic 
idea could either be proved or disproved with certainty because the method 
is based on the assumption that the airframe closely resembles a second 
order system. Objectives of the program were to determine: 

1) Will the iterative computing method converge to correct 
values of / and cOn ^ 

2) Conditions for convergence, 

3) Should both forms of the OJn equation be used, and if so 
V7hen is a particular form applicable? 

4) Effect of system damping ratio on convergence* 

5) Effect of sample spacing on convergence. 

The second order system used in the following programs was of the form: 

©^(DAMPiNG FACTOR)HI) 8-iSGS^O 

where the natural frequency == 6 radians per second was approximately the 
median of the airframe short period of this study, and the damping 
factor could be varied to produce damping ratios over the range from 0«1 
to i.O thus approximating the variation of the short-period mode damping 
ratio. Outputs were obtained by releasing the system from an initial dis¬ 
placement, 8 ^ at time = 0 and allowing it to oscillate freely. Samples of 
• • • 

0 , 0 , 6 c 0 v;ere periodically taken and stored for use in the / calculations. 
The dynamic simulation was used in four programs which were performed to 
test the damping ratio calculator and the influence of certain variables 
on its operation. The programs and results are discussed in the following 
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pages and sample listings of each are included in appendix VI, 


Program 1, Calculation of (A)r\ Using Known Values of ^ , 
and Calculation of ^ Using Known Value of • 


As a check on the accuracy of the dynamic simulationj> the samples of 
output were used in a series of non-iterative calculations in which U}f) 
was deteninined using the known value of^^and^ was determined using the 
known value of U)f) . 

It is recalled that the relationship for natural frequency, 



9 


e 


has two forms x^hich will be referred to for identification, as the posi¬ 
tive and negative forms. Both forms were used in t^ie calculations in order 


to determine under what conditions each V7as applicable. A third purpose 


of the program was to obtain a basis for comparing the accuracy of ^ and 
COf) computed by the iterative method. Nine runs were made in which damp¬ 
ing ratio was varied from 0.1 to 0.9 in O.l steps, and partial listings 
of the runs made at ^ =0.2 and 0.7 are presented in table IV, 

It is seen from table IV that the calculated values of ^ and cjp 

Q 

have approximately 0.5% accuracy except in regions where the ratios and 
0 

—g— become large. The errors in ^ and obtained with this series of 

runs are due to errors in the generated output values; 0, 0, & 0; and thus 

are indicative of the solution accuracy of the simulated second order system, 

and they also represent the highest accuracy which can be expected from the 

iterative method. It is also observed that the calculated values of 4 have 

very little variation xchereas the u/^ error increases in the regions where 

h (S 

0 is small, causing the ratios of —g— and to be large, and it is in 

t7 o 

these regions, before 9 changes sign, that the negative form of the u)n 
equation is applicable. Increasing damping ratio had the effect of widen¬ 
ing the regions in which OJp errors appear and in which the negative 
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Table IV Calculation of Natural Frequency Given Correct 
Damping Ratio, and Calculation of Damping Ratio Given Correct 
Natural Frequency for a Second Order System With CJn-6.0 Rad/ 
Sec and 3-0.2 and 0.7 

Sample Spacing 0.0125 Cycle. 


PLANT IDENTIFICATION ROUTINE - SECOND ORDER 

CONDITION= 2 DAMPING RATI0= .2000 


ACCELERATION 

V£ LOC I TY 

DISPLACEMENT 

OMEGAN1 

OMEGAN? 

OAMPRAT 

-3U.909819 

-. 434 735 

.997322 

6.004203 

-5.829842 

. 19.04 94 

-33.595382 

-.863053 

.989194 

6.004834 

-5.65584 1 

. 1946 19 

-32.133230 

-1.274004 

.975818 

6.005478 

-5.483248 

. 195985 

-30.535811 

- 1. 665819 

.957424 

6.006144 

-5.310185 

. 1966 72 

-28.815956 

-2. 036888 .934259 

6.006840 

-5. 134 753 

. 19 7089 

-26.986802 

-2.385762 

.906594 

6.007577 

-4.954950 

.197372 

-25.06U01 

-2. 711158 

.874713 

6.008367 

-4.768574 

.197578 

-23.054149 

-3.011961 

.838918 

6.009225 

-4.573107 

. 197736 

-20.977702 

-3. 287225 

.799521 

6.010 1 72 

-4.365574 

. 19 7863 

-18.845898 

-3. 536 1 73 

.756847 

6.011234 

-4. 14 2336 

. 19 7968 

-16.672187 

-3.758197 

.711229 

6.012447 

-3.898810 

. 1980 5 7 

-14.469854 

-3.952858 

.663006 

6.013863 

-3.629053 

. 198134 

-12.251956 

-4.119879 

.612522 

6.015557 

-3.325122 

. 198203 

-10.031248 ' 

-4.259145 

.560125 

6.017644 

-2.9760 75 

. 198264 

-7.820131 

-4.370700 

.506159 

6.020313 

-2.566302 

. 198 321 

-5.630584 

-4.454739 . 

.450972 

6.023889 

-2.072655 

. 198373 

-3.474119 

-4.511603 

.394904 

6.028996 

-1.459177 

. 198422 

-1.361722 

-4. 541775 

.338293 

6.036992 

-.666768 

. 1984 69 

.696184 

-4.545873 

.281469 

6.051495 

.408725 

.198514 

2.689788 

-4. 524639 

.224752 

6.086341 

1.966338 

. 198558 

4.609921 

-4.478935 

.168455 

6.272509 

4.362838 

. 198602 

6.448092 

-4.409733 

.112877 

9.794 1 77 

5.832562 • 

. 198645 

.. 8. 1965 13 

_-4.318107_ 

.058305 

23.690226 

5.934103 

.198688 

11.396610 

-4.072334 

-.046742 

5.972716 

-4 0.82 196’7 

. 198 7 78 

12.836410 

-3. ‘'2 0 762 

-.096718 

5.979735 

-22.194992 

. 198825 

14. 162/^31 
15.371551 

-3.751897 

-.144689 

5.984347 

-16.356640 

.198875 

-3.567184 

-,190449 

5.987638 

-13.479787' 

. 198929 

16.459617 

-3.368112 

-.233809 

5.990128 .. 

r11.752290 

.198986 

17.424440 

-3. 156207 

-.274598 

5.992097 

-10.589657 

.199049 

18.264289 

-2.933022 

-.312667 

5.993 7D8 

-9.745970- 

.199120 

18.978178 

-2.700125 

-.347884 

5.995064 

-9.099695 

.199200 

19.565846 

-2.459094 

-.380136 

5.996233 ’ 

-8.583825 

.199292 

20.027742 

-2.211504 

-.409334 

5.997261 

-8. 158339 

.199402 

20.364998 

-1.958920 

-.435403 

5.998182 

-7.797821 

.199535 

20.579it04 

-1.702891 

-.458292 

5.999021 

-7.485314 

.199704 

20.673379 

- 1.444938 

-.477967 

5.999795 

-7.209030 

.199925 

20.649939 

-1.186546 

-.494414 

6.000520 

-6.960482 

.200234 

20.512663 

-.929164 

-.507635 

6.001207 

-6.733358 

.200700 

20.265655 

-.674188 

-.517653 

6..00 18 66 

-6.522824 

.201496 

19.913508 

-. 422961 

-.524506 

6.002506 

-6.325066 

.203192 

19.461264 

-. 176 767 

-.528240 

6.003134. 

-6.136986 

.209474 

18.914369 

18.278638 

.063176 

-.52895 1 

6.003 758;- 

-5.955983 

. 168654 

.295721 

-.526699 

6.004383 

-5.779798 

. 192338 

17.560206 

.519797 

-.521593 

6.005016 

-5.606394 

. 195 131 

16.765485 

. 734408 

-.513744 

6.005666 

-5.433857 

. 196223 

15.901123 

.938644 

-.503276 

6.006340 

-5.260313 

. 1968 10 
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Table IV 


Continued 


ACCELERATION 

VELOCITY 

DISPLACEMENT 

OMEGANl 

_0.MEGAN2 

0AMPRAT 

14.973961 

1. 131675 

-.490325 

6.UU7U46 

-5.083842 

.19 7 180 

1 3.99U9d4 

1.312760 

-.475034 

6.007797 

-4.902394 

. 197437 

12.959284 

1.481246 

-.457558 

6.008604 

-4.713691 

. 197627 

1 1.8860 12 

1.636569 

-.438058 

6.009486 

-4.515100 

.197775 

10.778.344 

1.778253 

-.416701 

6.010463 

-4.303480 

.197895 

9.643431 

1. 9059 14 

-.393660 

6.01 1564 

-4.074955 

.197795 

8.488369 

2.019255 

-.369 1 13 

6.012829 

-3.824603 

.198080 

7.320 155 

2. 1 18068 

-.343239 

6.014315 

-3.545987 

.198155 

6.145654 

2.202232 

-.316222 

6.016107 

-3.230430 

. 19822 1 

4.971565 

2. 27 1 711 

-.288245 

6.018338 

-2.865862 

.193281 

3.804387 

2.326550 

-.259490 

6.021224 

-2..434886 

.198336 

2.650.39 1 

2.366876 

-.230141 

6.025154 

-1.911379 

. 198337 

1.515589 

2-392890 

-.200378 

6.030897 

- 1.254149 

.198436 

.405711 

2.404869 

-. 170378 

6.040203 

-.394233 

.193482 

-.673819 

2.403 160 

-. 140314 

6.058135 

. 792692 

.198527 

-1.717905 

2.388172 

-.110354 

6. 107544 

2.548848 

.198571 

-2.721799 

2.360380 

-.080663 

6.565521 

5. 139424 

. 1986 14 

-3.681113 

2.320313 

-.051396 

12.176150 

5.882218 

. 198657 

-4.591836 

2.268554 

-.022704 

34.023796 

5.9444 1 5 

. 198701 

-5.450344 

2.205735 

.005272 

5.964874 

-173.323327 

. 198745 

-6.253409 

2. 132527 

.032396 

5.975077 

-32.305450 

. 19879 1 

-6.998207 

2. 049643 .058545 

5.981237 

-19.985186 

.198839 

-7.682317 

1.957825 

.083600 

5.985395 

-15.352948 

. 198890 

-8.303729 

1.857846 

.107456 

5.988418 

-12.904145 

.198945 

-8.860837 

1. 750500 

.130016 

5.990736 

-11.376240 

.199004 

-9.352439 

1.636598 

.151191 

5.992589 

-10.322459 

.199069 

-9.777735 

1.516966 

.170907 

5.994 1 18 

-9.544512 

.199142 

-10. 136312 

1.392433 

.189095 

5.995415 

-8.940881 

.199225 

-10.428141 

1.263836 

.205 70 1 

5.996540 

-8.454162 

.199322 

-10.653565 

1. 132007 

.220678 

5.997534 

-8.049409 

.199437 

-10.813283 

.997771 

.233991 

5.998430 

-7.704088 

.1995 79 

-10.908339 

.861944 4 

.245615 

5.999248 

-7.402979 

. 199761 

-10.940106 

.725327 1 

.255536 

6.000007 

-7.135387 

.200003 

-10.910265 

. .588699 

.263748 

6.000720 

-6.893539• 

.200347 

-10.820793 

.452819 

.270257 

6.001398 

-6.67160.3 

.200881 

-10.673937 

.318418 

.275075 

6.002051 

-6.465079 

.201841 

-10.472195 

. 186 199 

.278226 

6.002687 

-6.270382 

.204106 

-10.218301 

.056831 

.279742 

6.003313 

-6.084575 

.216423 

-9.915192 

-.069053 

.-279662 

6.003936 

-5.905169' 

. 184 186 

-9.565996 

-.190857 

.278033 

6.004563 

•• '.-5.729981 

.193503 

-9. 174003 

-.308025 

.274909 

6.005280 

-5.557017 . 

.195530 

-8.742644 

-.420043 

.270353 

6.005856 

-5.384384- 

.196424 

-8.275467 

-.526441 

.264432 

6.006538 , 

. -5.210202 

. 196932 

-7.776116 

-. 626795 

.257218 

6.007256 

-5.032525 

.197262 

-7.248304 

-.720726 

.248789 

6.008021 

-4.849246 

.197497 

-6.695796 

-.807900 

.239228 

6.008848 

-4.658001 

.197673 

-6.122380 

-.888033 

.228621 

6.009755 

-4.456031 

. 1978 12 

-5.531850 

-.960888 

.217057 

6.010763 

-4.240008 

. 197925 

-4.927986 

- 1.026274 

.204630 

6.011907 

-4.005796 

.198020 

-4.314526 

-1.084048 

. 191432 

6.013228 

-3.748096 

.198102 

-3.695 155 

-1. 1341 12 

. 177561 

6.014792 

-3.459917 

.198174 

-3.073483 

-1. 176417 

.163112 

6.016693 

-3. 131756 

. 198239 

-2.453023 

- 1.210955 

.148182 

6.019084. 

-2.750264 

.198297 

-1.837183 

-1.237762 

.132870 

6.0222 187 

-2.295984 

.198351 

-1.229240 

-1.256917 

.117270 

6.026563 

-1.739318 

.198402 

-.632336 

-1.268539 

.101478 

6.033073 

- 1.032846 

. 198449 
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Table IV 
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C0NDITI0N= 7 DAMPING RAT 10= .7000 


ACCELERATION 
-32.518d30 
-29.005397 
-25.687L07 
-22.563123 
- 19.6299L3 
-16.884500 
-14.322759 
-11.940102 
-9.731413 
-7.691155 
-5.813440 
-4.092098 
-2.520738 
- 1.092801 
.198384 
1.359554 
2.397471 
3.318882 
4.130484 
4.838S95 
5.450624 
5.972048 
6.409392 
6.768709 
7.055865 
7.276529 
7.436159 
7.539995 
7.593052 
7.600117 
7.565747 
7.494262 
7.389753 

_ 7 oc: A ■> 70 _. 

* • j \j t r — 

7.096868 

6.915523 

6.715228, 

6.498949 

6.269440 

6.029253 

5. 780741 

5.526069 

5.267215 

5.005985 

4.744016 

4.482787 

4.223622 

3.967704 

3.716082 

3.469674 

3.229280 

2.995589 

2.769182 

2.550546 

2.340076 

2. 138085 

1.944809 


VELOCI TY 
-.420371 
-.804694 

- 1. 1 46321 

- 1. 447687 
-1.711196 

- 1.939218 
-2. 1 34075 
-2. 298034 
-2.433303 
-2.542022 
-2. 626260 
-2.688010 
-2.729188 
-2.751627 
-2. 7570 78 
-2.747209 
-2.723603 
-2. 687758 
-2. 64 1 089 
-2.584926 
-2.520520 
-2. 449037 
-2.371569 
-2.289128 
-2. 202652 
-2. 113008 
-2. 020993 
-1.927337 
-1.832705 
-1. 737703 
-1. 642875 

- 1.548 7 14 
-1.455656 

,,, 1 ? X ii n o ) 

I . W-» V.y / I 

- 1. 274361 

- 1. 186 762 
-1. 101552 
-1.018948 

-.939133 
-. 862256 
-.788436 
-.717763 
-. 650302 
-.586093 
-.525155 
-.467489 
-.413077 
-.361885 
-.313867 
-.268962 
-.227100 
-.188202 
-.152180 
-.118940 
-.088382 
-.060402 
-.034893 


DISPLACfcM 
.997379 
.989677 
.977440 
.961187 
.941405 
.918555 
.893063 
.865332 
.835732 
.8046 10 
.772284 
.739047 
.705169 
.670895 
.636449 
.602032 
.567826 
.533993 
.500677 
.468006 
.436089 
.405022 
.374888 
.345754 
.317676 
.290701 
.264861 
.240183 
.216662 
. 194366 
. 173238 
. 153292 
. 134516 

_ _ 1 1 A H O li 

• i t KJ / ^ 

.100406 
.085026 
.070727 
.057477 
.045242 
.033986 
.023672 
.014262 
.005715 
-.002009 
-.00895 1 
-.015152 
-.020652 
-.025492 
-.029712 
-.033352 
-.036449 
-.039042 
-.041 166 
-.042858 
-.044151 
-.045078 
-.04567 1 


NT OMEGANl 
6.012666 
6.012680 
6.012694 
6.012708 
6.01272? 
6.012737 
6.012753 
6.012769 
6.012786 
6.012804 
6.012823 
6.012843 
6.012865 
6.012890 
6.012917 
6.012947 
6.0 1298 1 
6.013020 
6.013065 
6.013120 
6.013186 
6.013269 
6.013377 
6.013524 
6.013737 
6.014078 
6.0147 1 1 
6.016315 
6.028037 
6.510341 
7.266745 
8.133335 
9.138630 

in zoi; c; AO 
• ^ ^ “ 

1 1 . 75 70 73 

13.528653 

15.792533 

18.807121 

23.049256 

29.507034 

40.6 162-90 

64.445604 

153.293882 

6.012401 

6.012427 

6.0-124 50 

6.0 124 72 

6.012492 

6.0125 10 

6.012527 

6.012544 

6.012560 

6.012575^.- 

6.012590' 

6.012604 

6.012618 

6.012632 


0MEGAN2 
-5.422600 
-4.874358 
-4.370803 
-3.904104 
-3.467937 
-3.057109 
-2.667295 
-2.294832 
-1.936570 

- 1.589 752 

- 1.251924 
-.92086 1 
-.59450 1 
-.270896 

.051839 
.375569 
.702179 
1.033626 
1.571978 
1.7J9973 
2.070580 
2.452076 
2.843 1 34 
3.255442 
3.693355 
4. 1620 70 
4.667840 
5.217939 
5.813242 
6.006137 
6.0099 10 
6.010925 
6.011397 

A n 1 1 A 7 o • 

4 / • V./ I I \ / I 

6.0 1 185-2 
6.0 1 1979 
6.012075 
6.012150- 
^ 6.012211 
^ 6.012261 
6.012304., 
6.0 1 2341 
6.012373 
■414.409029 
-88.150044 
-49.208001 
-34.015191 
-25.886909 
-20.801486 
-1 7.302693 
-14.735417 
-12.761324 
-11. 187999 
-9.897905 
-8.815159 
-7.888553 
-7.082255 


DAMPRAT 

.671394 

.685868 

.690646 

.693036 

.694 4 78 

.695446 

.696 145 

.696676 

.697096 

.69 74 37 

.697721 

.697963 

.698173 

.698358 

.698522 

.698670 

.698805 

.698729 

.699045 

.699153 

.699254 

.699 3 5 1 

.699444 

.699533 

.6996 19 

.699703 

.699786 

.699867 

.699948 

.700029 

.700 1 10 

.70019? 

.700275 

/nn ^ A 1 

• I \J KJ yy I 

. 700448 
. 700538 
.700633 
. 700 73 1 
.700836 
. 700947 
.701066 
.701194 
. 70 1335 
.701489 
.701662 
.701856 
. 7020 79 
.702338 
.702645 
.703017 
.703479 
. 704074 
.70487? 
.706010 
. 70 7 7 73 
.710899 
.718024 


( 
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Table IV 
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\ 


CCELhRATION 

VELOCITY DISPLACEMENT 

' OMEGANl 

0MEGAN2 

DAMPRAT 

1.760in2 

-.011 745 

-.045960 

6.0 12645 

-6.3704 15 

.75 10 n 

1.5B4997 

.009154 

-.045974 

6.012659 

-5.73389 1 

.637835 

1.418607 

.027917 

-.045740 

6.012673 

-5. 158183 

.680<-.93 

1.261231 

.044657 

-.045285 

6.012687 

-4.632083 

. 6880 O 2 

1.112814 

. 059486 

-.044632 

6.012701 

-4. 146765 

.691943 

.973255 

.072514 

-.043805 

6.012715 

-3.695161 

.693794 

. .342418 

.083853 

-.042826 

6.012729 

-3.271520 . 

.694977 

.720133 

.093611 

-.041715 

6..012745 

-2.871087 

.695802 

.606202 

.101892 

-.040492 

6.012760 

-2.489869 

.6964 12 

.500400 

. 108799 

-.039 1 74 

6.012777 

-2. 124461 

.696885 

.402483 

.114434 

-.037777 

6.012794 

-1.771913 

.697265 

.312188 

.118893 

-.036318 

6.012812 

-1.429621 

.697577 

.229238 

.122270 

-.034809 

6.012832 

-1.095246 

.697840 

.153343 

. 124654 

-.033265 

6.012853 

-.766647 

.698066 

.084206 

. 126132 -.031697 

. 126786*"^^ -.0301 15 

6.012877 

-.44 1823 

.698263 

.021523 

6.012902 

-. 1 18860 

.698437 

-.035015 

.126695 

-.028530 

6.012930 

.204 1 1 1 

.698594 

-.085719 

.125935 

-.026951 

6.012962 

.528959 

.698735 

-.130900 

.124575 

-.025384 

6.012998 

.857599 

.698865 

-. 1 70867 

.122684 

-.023838 

6.013040 

1.192032 

.698985 

-.205927 

.120324 

-.022319 

6.013090 

1.534400 

.699096 

-.236382 

.117555 

-.020832 

6.013149 

1.887041 

.699201 

-.262529 

.114432 

-.019382 

6.013222 

2.252559 

.69930 1 

-.284657 

.111009 

-.017972 

6.0133 16 

2.633904 

.699395 

-.303050 

.107332 , 

r -.016608 

6.013440 

3.034477 

.699486 

-.317981 

.103447 

-.015290 

6.013614 

3.458259 

.6995 74 

-.329718 

.099395 

-.014022 

6.013876 

3.9099 75 

.699659 

-.338515 

.095216 

-.012806 

6.014320 

4.395298 

.699742 

-.344620 

.090944 

-.011642 

6.015248 

4.921011 

.699824 

-.348268 

.086611 

-.010532 

6.018406 

5.494219 

.699906 

-.349685 

.082246 

-.009477 

6.158836 

5.991120 

.699786 

-.349086 

.077877 

-.008476 

6.854105 

6.008693 

.700067 

-.346677 

• . 0 73527 

-.007530 

7.659854 

6.010503 • 

.700149 

-.342650 

.069217 

-.006638 

8.587396 

6.01 1 184 

.700231 

-.337189 

.064966 

-.005799 

9.671900 

6.011543 

.700315 

-.330467 

.060792 

-.005013 

10.964613 

6.01 1765 

.700401 

-.322.644 

.056709 

-.004279 

12.541722 

6.011917 , 

.700490 

-.313874 

.052730 

-.003595 

14.52 1378 

6.012028 

. 700582 

-.304298 

.048866 

-.002960 

17.097171 

6.0121 13 

. 700679 

-.294047 

.045126 

-.002373 

20.609855 

6.012180 

.700780 

-.283244 

.041517 

-.001832 

25.7201 14 

6.012236' 

. 700887 

-.272002 

.038046 

-.001335 

33.899489 

- 6.012282 

.701002 

-.260425 

.034718 

-.000080 

49.225428 

6.012322 

.701 125 

-.248609 

.031537 

—.000466 

88.7339 17 

6.012356 

.701259 

-.236640 

.028504 

-.00009 1 

432.99657-9 

6.012387 

.701406 

-.224599 

.025621 

.000247 

6.0'T24 14 

-151.100646 

.701568 

-.212556 

.022889 

.000550 

6.012438 

-64.247932 

.701750 

-.200577 

.020307 

.000820 

6.012461 

-40.679637 

.701957 

-. 188719 

.017874 

.001059 

6.012481 

-29.65 1905 

.702196 

-.177035 

.015588 

.001268 

6.012500 

-23.2297 14 

.702476 

-.165568 

.013447 

.001449 

6.012518 

-19.006207 

.7028 1 1 

-.154359 

.011448 

.00 1604 

6.012535 

-16.002548 

.703221 

-. 143443 

.009587 

.001736 

6.012552 

-13.745592 

.703739 

-. 132848 

.007860 

.001845 

6.012567 

-1 1.978633 

.704418 

-.122599 

.006264 

.001933 

6.012582 

-10.550309 

.705354 

-.112717 

.004794 

.002002 

6.012596 

-9.365561 

.706735 
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equation is applicable* It >ouM be expected from the results of this 
series of runs that the iterative method of determining ^ and ojp niight 

experience difficulty in converging whan G is small and that the region of 

non convergence will become wider as system damping ratio is increasedo 

Program 2. Iterative Calculation of ^ and Using 
Both Forms of the Equation. 

The purpose of this run was to test the iterative damping ratio com¬ 
puter described in section V. The previously described second order system, 

• ft % 

with €* & ^=6>.2. , was used to supply output samples of 0, 0, and 0 

which were stored, two sets at a time, and used as inputs to the comput¬ 
er. Sample spacing in this run was 0.05 cycle and each calculation 

was iterated 20 times using each of the two forms of the OJn equation. At 

the termination of each iterating process, the results of each of the 20 
iterations was printed out for study, the oldest set of samples of 0, 0, & 

0, discarded, a new set generated and stored, and the calculation cycle 

repeated using the two new sample sets. A sample series of the printouts 
from the computer is listed in Table V where the results of one con¬ 

vergent and one divergent solution are shovm. 

The primary result of this program was the demonstration that the itera¬ 
tive method would converge, and it did so for* 51% of the sample sets - 
0.2). Accuracy obtained w^as of the same order as that obtained in Program i?l, 
approximately 0.05%, and in regions of small values of 0 the accuracy of the 
iterative method was much greater. As an example, from the iterative 
method was 6.00359 as compared to ss 6.086341 obtained from program 
when the identical sample set v/as used with each method in a region of 
small values of G. It v;as found that a necessary but not sufficient condi¬ 
tion for convergence of the iterative method v/as that 0 and 9 be of opposite 






As is shown in Tnble \\ another interesting property of the iterative 
method is that when convergence occurs, it occurs simultaneously with both 
forms of the equation, and the iterative series using the negative form 
of the cJfj equation converges to the negative values of ^ and arrived 
at using the positive equation. Also, it was found that the use of 
the negative equation would not be required in future programs because 
the iterations performed when using the positive form converged in areas 

where the results of Program #1 indicated that the negative ^ equation was 

0 

applicable. It had thus far been established that, at ^ ^ 0.2, the itera¬ 
tive computations would converge for a length of time corresponding to approxi¬ 
mately one quarter of the system period at a rate of two times per system 
cycle. It was felt that this would suffice for the proposed adaptive 
system, however, it was considered desirable to determine the effect of 
different values of damping ratio on convergence. 

Program 3. Iterative Calculation of / at 9 

Values of System Damping Ratio. 

Nine runs were made in this program, with system damping ratio varying 
from 0.1 to 0.9 in 0.1 steps, and a sampling interval ^ 0.05 period. Only 
the positive form of the equation v/as used in the Iterative process and 
the solution was tested for divergence after each iteration to prevent un¬ 
necessary calculations. Table VI is a partial listing of the runs made with 
/ = 0.2 and 0.7. It is seen that the region of divergence is appreciably 

wider with / =0.7 than with 0.2. The percentage of convergent iterations 

to total iterations is tabulated below as a function of damping ratio for 
each of the nine conditions; 
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RStlo/uslnrBotrForarorthe ^ Frequency and Bampins 

System Klth®<«,. 6.0 Rad/Sec L"'’/^T2 ’ “ 

Sample Spacing = 0.05 Cycle 


TIKE 

0,2ccccr) 




\ 


-Fo 


or/T) 


ACCELERATION VELOCITY DISPLACEMENT 
-5o63v577 -4o45473G 0o45!'i97T 


omegan 

7o629995 
-- 6,818375 
6,413184 
6,21ur'9] 
6, U’789G- 
6,056324 
6,^T3>)2o7 
6,017089 
6«')lu425 
6,307049 
6,0-”-5342 
6,0C4477 
6,004040 
6,003819 
6, 003707 
6, }•: 3650 
6,003621 
6,0(5607 
6onr-5599 
6,0r-559£ 
OMEGAN 
-2,860373 
-4,239639 
-5.'‘7871 1 

-5,528898 

-5,761842 

-5,880874- 

-5.941390 

-5,972086 

-5,987659 

-5,995516 

-5,999504 

-6,'51522 

-6,06-2544 

— 6 o 03 3 6 2 
- 6,:'':3323 

- 6 , C”0 3 4 5 6 
-6,-3C3523 
-6,003557 ■ 
-6,003574 
-6.003583 


OAMPR 

0,303380 
0,252439 
0,226073 
3,212571 
0,205695 
0,202204 
0,23:'434 
0,199537 
1990 83 
0,198854 
C,198737 
7,198678 
0, 198648 
‘•o 198633 
0,198626 
0,198622 
0,198620 
Oo198619 
0,198618 
0, 198616 
DAM PR 

0,r*761 £8 
-0.065534 
-0,132633 
-0,165552 
-3,181965 
-3.19 021./ 
r ■>o r94367 
-0,196468 
-3<rl9753:. 
-0,198067 
rOo198339 
-0,198477 
-0,198547 
-0,198582 
-Oo 1 986 CO 
-0, 198609 
-:u; 1986 13 
-0,198616 
-0,198617 
-0,198617 


9 


/ f 
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Table V - Continued 


f 

TIf'E ACCELERATION VELOCITY DISPLACEMENT 

0.399999_1 8,978094__ -2.7001.15_-0.347882'.„_ 


OMEGAN 

DAMPR 

4.205901 

0.564624 

3.300829 

0.852037 

2.354981 

1.340582 

1.558809 

2.154067 

0.988132 

3.492864 

0.613964 

5.684419 

0.378345 

9-264260 

0.232391 

15.107418 

0.142551 

24.643762 

0.087360 

40.222131 

0-053459 

65.734811 

0.032669 

107.57C982 

0.019852 

177.027143 

0.011444 

307.084564 

_ 0.005859 

599.775604 

n.002441 

1439.462204 

-0.002930 

-1199.551773 

44936.139648 

-2894.778473 

108440.639648 

-6985.727661 

261690.726562 ' 

-16858.072266 

OMEGAN 

DAMPR 

-19.190211 

1.G531C0 

-41.399521 

2.582062 

-97.553423 

6.248351 

-234.412361 

15.085813 

-565.269005 

36.408317 

-1363.942886 

87.862394 

- -3291.417694 

212.031506 . 

-7942.874573 

511.678204 ^ 

-19167.864990 

1234.790436 

-46256.204590 

2979.816833 

-111626.239258 

7190.943481 

-269378.281250- 

17353.303223 

-650068.140625, 

41877.279297 

-1568755.15625C 

1C1058.9414Q6 

-3785745.781250 

243877.103516 

-9135824.37500C 

588528.257812 

-22&4672C.500000 • 

1420246.000000 


i4S7360«906SS5»_„. 

-128391576.000000 

8270963.625000 

_ _-77459164.000^00_ 

19959625,000000 _ 
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Table VT 


- er .:er«t o ^o^ e Tterations as a Function of Syster 

Tar.piu;^ ^.T^io 


Damping Ratio 

O.l 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

7o Convergence 

52% 

517o 

51% 

49% 

48% 

41% 

34% 

27% 

19% 


The percentage of convergent Iterations drops off sharply for values of 
exceeding 0.5 but it is not desired to increase the airframe damping ratio 
above 0.7 and thus the minimum obtainable ratio of convergent to total itera¬ 
tions should never be less than approximately 1/3. 

Program 4. Influence of Sample Spacing on the Percentage 
of Convergent Iterations. 

The final program was run to determine what influence^ if any^ the width 
of sample spacing has on the ability of the iterative calculations to converge. 
All runs were made with / “0.2, and six values of sample spacing ranging 
from O.Ol cycle to 0.5 cycle were investigated. Results of the program are 
tabulated below. 

Table VIIT - Percent of Convergent Iterative Calculations as a Function of 
Sample Spacing / =0.2 


Spacing 

0.01 cycle 

0.02 

0.05 

0,1 

0.2 

0,5 

% Convergence 

56% 

53% 

51% 

50% 

48% 

46% 


It is seen that narrow sample spacing favors convergence^ however it can be 
varied over wide limits vzith little effect on convergence. It was arbitrar¬ 
ily decided to set sarple spacing at 0.025 seconds for the three degiee of 
freedom airframe simulation; the primary consideration being conservation of 
computer time. This sample spacing corresponded to approximately 50 samples 
per short period cycle, depending on flight condition. 
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Table VII Iterative Calculation of Natural Frequency and Damping 
Ratio of a Second Order System With Rad/Sec and/= 0.2 ^ 0.7 

Sample Spacing = 0.05 Cycle 


Condition 2 

Damping Ratio « 

0.2 


Tlf^E 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.2C0000 

-5.630577 

-4.454730 

0.450971 

Of'EGA-N 

OAMPRAT _ 



7.629993 

0.303380 



6.818375- 

_0.252439 .. 



6.413184 

0.226073 



6.210091 

0.212571 



6.107890 

0.205695 



6.056324 

0.2C2204 



6.030267 

0.200434 



6.017089- 

0.199537 



6.010423 

0.199083 



6.007049 

0.198854 _ 



6.005342 

0.198737 



6.004477 

0.198678 



6.004040 

0.198648 



6.003819 

0.198633 -- 



6.003707 

0.198626 



6.003650 

0.198622 



6.003621 

0.198620 



6.003607 

0.198619 



6.003599 

0.198618 




6.003595 0.198618 


TIME ' 

ACCELERATION 

VELOCITY DISPLACEMENT 

0.250000 

2.689775 

-4.524627 0.224752 • 

OMEGA-N 

DAHPRAT 

' 

9.251958^ 

0.261913 


6.966586 

0.215692 

*• N. 

6.256753 

0.202902 


6.066626 

0.199669 


6.019025 

0.198874 

•• 

6.007354 

0.198680 


6.004508 

0.198633_ 


6.003815 

0.198622 


6.003646 

0.198619_ 


6.003605 

0.198618 


6.003595 

0.198618, 


6.003592 

0.198618 •• 


6.003592 

0.198618 


6.003592 

0.198618 


6.003592 

0.198618 


6.003592 

0.198618 


6.003592 

0. 198618 


6.003592 

' 0.198618 
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Table VII- Continued 


Damping Ratio - 0.2 




i 



T 


6.003592 0.198618 

6.003592_ 0.198618 


TIME 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.299999 

_ 9.848088 

-4.205210 

0.005013 

OMEGA-N 

DAMPRAT 



___15.324385 . 

_ 0.085544 _ 



-7.279872 

-0.165185 



.-4.234596 

-0.279041 _ _ 



-1.191589 

-0.983381 




-SOLUTION DIVERGING 


TIME 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.349999_.. 

15.371489 

-3.567171 

-0.190448 

OMEGA-N 

DAMPRAT 



SOLUTION DIVERGING 







TIME 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.39999-9 

_L8.-9-7809Jt_ 

-2.7QQ115 

-0.347882 


OMEGA-N DAMPRAT 

_4.205901__0.564624 

3.300829 0.852037 

-SOLUTION DIVERGING _ 


TIME ACCELERATION VELOCITY DISPLACEMENT 

_0.449999_20.579304_^1.702884 __-0.458290_ 

OMEGA-N DAMPRAT 

_4.907358 _0.570963__ 

4.181908 • 0.882182 

-SOLUTION--OI-VERGING- 


TIME ACCELERATION VELOCITY DISPLACEMENT 

_0.499998_20.265548 — -0.674185 _-0.517650_ 

OMEGA-N DAMPRAT 

__5.37J 64 3_0 .-130322_ 

SOLUTION DIVERGING 


_ TIME _ACCELERATION_^VELOCITY_DISPLACEMEN1_ 

0.549998 18.278530 0.295719 -0.526696 

_QMEGA-N D AMPRAT_ * _ 

5.757618 -0.240352 

SOLUTION-DIVERGING___ 
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Table VII- Continued 


Damping Ratio 


0.2 


TINE 

ACCELERATION 

VELOCITY DISPLA-CEMENT 

0.599997 

14.973863 _ _ 

1.131666 _-0.490321_ 

QNEGA-N 

6.119884 

DAMPRAT 

0.244751 


6.030041 

0.209181 


6.009638 

0.201036 


6.004975 

0.199172 


6.003908 

0.198745 


6.003664 

0.198647 


6.003608 

0.198625 


6.003595 

0.198619 


6.003593 

0.198618 . _ 


6.003592 

0.198618 


6.003592 

0.198618 


6.003592 

0.198618 

. 

6.003592 

0.198618 


6.C03592 

0.198618 


6.003592 

0.198618 


6.003592 

0.198618 


6.003592 

0.198618 


6.003592 

0.198618 


6.003592 

0.198618 


-6.003592 _ 

0.198618 



TINE 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

0.649996 

10.778271_ 

1.778240 

-0.416698- _ 

ONEGA-N 

DAMPRAT 



6.525987. 

0.300233 



6.262415 

0.249806 



6.132749 

0.224382 _ 



6.068288 

0.211580 



6.036061 

0.205138_ 



6.019904 

0.20189.7 


... - 

6.011790 - 

0.200267_ 



6.007714 

0.199447 

*• 


-6.005664 

0.199035 



6.004634 

0.198828 



.6.004116 

0.198723., 



6.003855 

0.198671 


. * 

6.003724. 

0-198645 



6.003658 

C.198631 



6.003625 

- 0.198625_ 



6.003608 

0.198621 



6.003600 

0.198620 



6.003596 

0.198619 



6.003594 

0-198618 



6.003593 

0.198618 













































I 


Table VH- Continued 

Condition 7 Damping Ratio * 0.7 


TINE 

^ACCELERATION VELOCITY 
_ 1.359547 -2.747204 

DISPLACEMENT 

_Q.20000.Q . 

0.602031 

ONEGA-N 

DAMPRAT 


„4.221364^ 

0.521158 


4.917109 

0.589098 


5.Z25D53 _ 

0.629943 


5.575466 

0.655293 


_5.732646 

0.671299 


5.832533 

C.681505 


5.896477 

0.688051 


5.937594 

0.692266 


5.964107 

0.694986 


5.981233 

C.696743 


_5.992308 

C.697880 ' 


5.999476 

0.698616 


6.004117 

0.699093 


6.0C7122 

G.699402 


6.009069 

C.6996C2.. 


6.010331 

C.699731 


6.011 148 

0.699815 


6.011677 

0.699870 


6.012021 

C.699905 


6.012243 

C.699928 



TINE 

ACCELERATION VELOCITY 

DISPLACEMENT 

0.25OCOO 

4.838883 -2.584920 

0.468005 

ONEGA-N 

DANPRAT 


2.861347 

0.586139 


4..8.ej„^9J 

r.633940 


5.364674 

0.660114 


5.622871 

C.675476 


5.773553 

0.684772 


5.864458 

0.690488 _ 


5.920257 

0.694035 


_5.954.843 

C.696248 . 


5.976406 

0.697633 


5.989896 

(5.698501' . 


5.998354 

0.699046’ 


6.003663 

0.699389 


6.007000 

0.699605 


6.0C9097 

0.699740 . 


6.010416 

0.699825 

• 

6.011245 

0.699879 


6.011767 

0.699913 


6.012C96 

0.699934 
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Table VII- Continued 


Damping Ratio 


0.7 


6.012302 

6.012432 

u,699947 
699956 

— 


TINE 

ACCELEPnATION 

VELOCITY 

DISPLACEMENT 

0.299999 6.768687 

ONEGA-N DAMPRAT 

-3.249:22 -•J.7'0412 

_ -1.717651 -0.990456_ 

SCLUTICN DIVERGING 

-'’.289121 

» 

0.345753 

TINE 

0.349999 

ONEGA-N 

-9.915053 

-2.3C6982 

ACCELERATION 
7.539968 
__ DAMPRAT .. 
-0.815083 
-C.991656 

VELOCITY 
-T.927330 

DISPLACEMENT 

0.240182 

SOLUTION DIVERGING 



TINE 

0.399999 
. .. ONFGA-N 

ACCELERATION 
7.494233 
.C.4MPRAT_ 

VELOCITY 

-1.548708 

DISPLACEMENT 

153291 

SOLUTION DIVERGING 

_ TINE __ 

ACCELERATION 

velocity 

DISPLACEMENT 

0.449999 
_ ONFGA-W . 

6.915493 

OANPRAT 

-•'. 186757 

0.085026 

SOLUTION DIVERGING 

TINE 

ACCELERAT ION__ VELOCITY 

DISPLACEMENT 

0.499998 
nNFDA-N 

6.029223 
nAMPRAT: 

-P.862252 

0.0339 86 

-44.580629 -0.957002 

SOLUTION DIVERGING^ 







TINE 

. _.C.549998 

ACCELERATION 
5.Of 5958 

VELOC ITY 
-P.586093 

DISPLACEMENT 
• ' -0,0020C9 

CNEGA-N 

-64.266227 

DAMPRAT 

. G.043701_ 



SOLUTION DIVE 

RGING 



_TINE 

. ACC.EL£8_AT ION 

VELOCITY 

DISPLACEMENT 

0.599997 
ONEGA-N 
IG.228298 
20.289338 

3.967681 

DAMPRAT 

0,175711 
444423 

-'',361 883 

-0.C25492 


SCLLTICN DIVERGING 


71 


















































I 


Table VII- Continued Daraplns Ratio » 0.7 > 


_TII'E 

0.649996 

OVEGA-N 

ACCELEBATION, 

2.995569 

CJVNPRAT 

0.158457 

-0.316348 

_VELOCiTY_ 

-''.18820C 

DISPLACEMENT 

-0.039041 

8.C28860 
_10.427491 



SOLUTION DIVE 

RG IN G 



TINE- 

ACCELERATION 

VELOCITY 

DISPLACEMENT 

1.:^ 0.699996 

2. 138C70 

.0 6040'2__„ 

-0.045078 

ONEGA-N 

DAKPRAT 



7.040965 

-C.113642 



SOLUTION DIVERGING 


TIME 

ACCELERATION VELOCITY 

DISPLACEMENT 

0.749995 

1.418596 027917 

-0.045740 

OMEGA-N _ 

DAMPRAT 


SOLUTION DIVERGING 


TIME 

ACCELERATION VELOCITY 

DISPLACEMENT 

0.799994 

0.84241 1 ''.083853 

-0.042826 

CMEGA-N 

DAMPRAT 

• 

5.818552 

0.622535 


5.961974 

0.679927 


5.999497 

0.694778 


6.CC9243 

0.698625 ... 


6.011769 

0.699621 


_ 6.012424 

0.69988C 


6.012594 

C.699946 


6.012637 

0.699964 


6.012649 

0.699966 


6.012652 

0.699969 


6.012653 

0.69997C 


6.012653 

u.699970 


6.012653 

C.699970 


_6.012653 

0.699970- 


6.012653 

0.699970 


6.012653 

... C.69997U._ - . . ... 


6.012653 

0.699970 


_ 6.012653 

0.699’97=) 


6.012653 

0.699970 
fl. 69-99.21) 
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VII “ APPLICATION Op THL OAMPTNG KiJfTO COfirUTEP TO A SELF ADAPTIVE PITCH 

DAMI IMG SVSTEM 

Investigation of the Ability of the Computer to Determine 

Short-Period Mode Characteristics 

It has thus far been established that the damping ratio computer opei- 
ates satisfactorily when supplied input samples generated by a second order 
system, and it has also been established that the short period mode, a 
second order oscillation, is the dominant mode of the airplane longitudinal 
dynamics. The next step in the program is then to investigate the opex'a- 
tion of the computer when sut plied samples generated by the three degree of 
freedom, fourth order airframe simulation. The simulation program is dis¬ 
cussed in detail in Appendix VI. 

Flight condition number seven^M =: 1.8 25,000 ft., having short-period 

mode damping ratio = 0.095 and period = 1.2 seconds, was chosen for the test. 
At the start of the computer run, the airframe x^as released from a non-trim 
pitch condition and allov/ed to oscillate. Periodic samples of pitch angle, 

0, and its first two derivatives v;ere taken and supplied to the com¬ 

puter as described in the proceeding section. At the beginning of each set 
of calculation an initial value of ^ = 0.4 x^as introduced into the 

equation and the iterations w^ere carried out 15 times. The iterations con¬ 
verged for 49 percent of the sample sets, x^hich is in good agreement with 
the previous second order system investigation, hex^ever the computed values 
of damping ratio ranged from 0.1448 to 0.024565 during a 2.5 second (flight 
time) run. Such a wide variation of computed ^ v/as inconsistent with the 
results of the preceeding study in which it was demonstrated that the longi¬ 
tudinal aircraft dynamics very closely approximate the dynamics of a second 
order system. A reviex^ of foregoing investigations x^as conducted to deter¬ 
mine if there x>7ere any significant differences between the second order 
study and the present application of the computer^ and the answer \^as 
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found in c xa^uln ‘ t j. .n «. •* atrriijr.? sbort-poriod pltrb angle-to-c levator 
deflect ion t rans Fer furtcf 1 

j9 ^ /r. (6^Te) _ 

dCS ^ C^nS +‘ 

It is seen that the 1/S tenri causes pitch angle output to exhibit a theo¬ 
retically infinite steady state error or^ if the alifrarne is perturbed froin 
its trim position at a given f^light conditionj^ the ensuring short period oscil 
lations will not decay to the original trim angle. Further^, it can be shown 
that the airframe will return to trim condition as a function of the long}- 
tudinal phugoid mode. The primary difference betv;een the ideal second order 
system and the short period dynamics is then the fact that the former system 
oscillates about 0 equal zero whereas the airframe oscillates about an aver¬ 
age pitch angle which is constantly varying as a function of the phugoid 
dynamics, and thus the pitch angle samples introduced into the computer were 
erroneous. The solution to this problem was seen in the pitch-rate-per-stabi 
lizer-deflection transfer function, 

Q _ _ (Si' Ts ) _ 

't (6^-h^^U)n6-h 0)n^) 

which exhibits zero steady state error, and which also shovjs that free air¬ 
frame short period oscillations oscillate about, and decay to, zero pitch 
rate . It was realized that pitch rate is the true second order variable of 
the airframe longitudinal dynamics and in the study to be outlined next, 
pitch rate and its first two derivatives 9 and G were utilized as inputs to 
the ^-uln computer. 

A second test of the computer using the three degree oc freedom 
flight simulation was performed, as stated above, and a partial listing of 
the results obtained when sampled values of 0, 0 0, were used as computer 

inputs is presented in Table IX. As before, the calculations converged for 
49X of the sample sets but, as can be seen, the calculated values of damp¬ 
ing ratio compare favorably with the known ^ 5 = 0.095. It can also be seen 
from examination of Table IX that, in most rases, the. calculations have not 
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quite settled ou^ to srcv«. iy \ i! 3S • c ^ and aft' 15 i ttrar ionr>^ and 
as a resultj it was decided to re[eat the itetatlve ca'calation cycle 20 
tiires in all future pcograniSo 

Application of the computer to the Adaptation Loop 

It was felt^ at this time^ that the computer was generating suf¬ 

ficiently accurate values of damping ratio to allow its use in the adapta¬ 
tion loopj how^ever one obvious discrepancy remained to be corrected before 
the loop could be closed* That iSj, what should be done regarding the 
calculations v^hich do not converge and thereby generate highly inaccurate 
values of ^ ? Before discussing the technique which \;as used to detect and 
account for divergent calculations^ possibly a better understanding may be 
obtained from a short review of the proposed mechanization of the adapta¬ 
tion loop. 

It x^ill be recalled that it was proposed to compare the output of the 
damping ratio computer with the desired value of short period damping ratio 
j the error between ^ - desired and ^- computer being used as the input 
signal to a rate servo which would vary pitch damper gain in the proper 
direction to null out the error, and at a rate proportional to the 
magnitude of the error* It is apparent that the erroneous values of ^ 
resulting from divergent calculations could not be allowed to reach the 
computer output, and it was found that this could be accomplished by 
simply including a routine which examined the resultant values of ^ 

at the termination of the series of iterative calculations for each sample 
set. After examination of a large number of past output records from both, 
the second order, and the airframe studies, it was determined that, in 
every case, when the calculations diverged the r^esultant calculated values 
of either ^ or or both were completely out of the range of variation 
of the real values of ^ and * For examples of divergent calculations, 
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Table IX Iteraoive Calculation of Airframe Short Period Damping 

Frequency Using 9, 6, & 0', Generated by a iLg- 
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Table IX 
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see Tables V ov VII, n'^so, t v;as Unovn with certainty, that for the air¬ 
craft considered in this study, damping ratio was restricted well within 
the range, from C\ to 1,0 and natural cx'equency within the range from 0. 
to 8,0 radians/second* Thus a divergence detection routine was programmed 
to use a simple logical process for determining whether or not the calculated 
values of both ^ and CJr) fell within the allowable ranges« If both values 
V7ere acceptable to the detection routine, the calculated value of ^ was 
set into the damping ratio comparison process, and if either ^ or or 
both were unacceptable to the detection routine, the last acceptable com¬ 
puted value of ^ was reset into the comparator. The values of 0)^ and ^ 
resulting from a divergent series were always orders of magnitude different 
from the allovrable ranges of uJf^ and ^ , and thus the allowable ranges 

could be altered considerably to suit the characteristics of other systems 
with no deterioration in the effectiveness of the divergence detection 
routine* As an example, if it were expected that the airframe might have 
short-period poles in the right half plane over some flight conditions, the 
alloxcable ^ range could be increased from ^ = - l.O to + 1*0. 

A variable v/hich has not yet been discussed, but v^hich was quite im¬ 
portant in the evaluation of the adaptation loop, is the gain of the rate 
servo (hereinafter called the adaptation servo) used to vary pitch damper 
gain* When ^ - error* signals are supplied to the adaptation servo which 
in turn causes variations in airframe damping ratio, the entire adaptation 
loop appears to resemble a second order servomechanism and should be amenable 
to conventional systems analysis techniques. Hox/ever analysis was compli¬ 
cated by the intermittent action of the damping ratio computer and thus an 
’’optimum*’ value of adaptation-seivc-gain, K^, was determined experimentally 
by performing a series of three degree of freedom computer runs with vary¬ 
ing values of Flight condition 7, M = 1.8 ''^.25,000 ft,, \ms again used 


78 


H': 







Figure 29. Time History of Airplane & Control Surface Motion 
and Pitch Damper Gain During Adaptation Process Vn’ith No Adapt¬ 
ation Servo Lock 

Initial = 0, Adaptation servo gain = 1. 
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for r.ll runs and a.IaplMtion-servo gain v/as vax^led fcoui K ^ 0.05 to 2.0 

s 

in six steps. 

Figure 29 is a typical time history of the run made v;ith K «= 1.0, 

s 

Airplane pitch angle, pitch rate^ control surface position^ and pitch damp¬ 
er gain are shovrn as functions of time^ and again the run V7as started by re¬ 
leasing the airfrarrje from a non-trim pitch position and allowing it to oscil 
late. The damping ratio of the basic airplane at this condition was 0.096 
and the steady state damper gain for correct adaptation^ from the root locus 
plot, was 0,344. Initial damper gain was set at zero at time ~ zero and it 

is seen that K, = 0.344 was achieved in 0,45 seconds, and that the airframe 
a ^ 

motion w^as heavily damped by the control action. However, it is also seen 

that damper gain did not level out at the desired value of K, = 0,344 but in 

d 

stead continued increasing In a somewhat random manner until at time ss 10 
seconds (not shovm) it reached the required value for instability of the 
actuator mode and a violent divergent oscillation ensued. In this condition 
the actuator motions comprised the dominant longitudinal mode, and it was 
interesting to note that the computer was still performing its duty in 

that it was now computing natural frequency and damping ratio of the actua¬ 
tor mode. The crossover frequency for this mode was predicted to be 34 
radians per second from the root locus plot, which was verified by the fact 
that computed damping ratio changed from positive to negative sign v;hen 
computed Ci)^ v;as approximately 34.5 rad/sec. But, the operation of the 
adaptation loop was obviously unacceptable with = 1.0, An improvement 
in performance was detected when loiv^er values of adaptatlon-servo-gain were 
Investigated, and with K, ~ 0.1 the damper gain did settle out reasonably 

a 

well, but approximately eight seconds v/ere required to complete the adapta¬ 
tion process when using the K^, Osburn has shevm in Ref. 127 that the MIT 
self adaptive system is capable of completing the adaptation process on a 
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siri.llr v nirfr/ .o. lurin;-^ lv..o v.yc1es or \ osc 11 lac wnich 'i^ovild 

come^>j“Oncl to ada’; tatlon tine - 2,4 secorrls for flight condition 7, 
presently Jiscussion. Thus, a reviev; van made of the adaptation 

nethod to detemiine if the process could be speeded ur without attend¬ 
ant drift in pitch damper gain. 

The inability of the adaptation loop to sec proper damper gain when 
operating at high adaptation-servo-galns was fount! to be caused by foTir 
interacting factors: 


The reader will recall that when the process failed to 

converge, the computer set the most recently comp^uted con¬ 
vergent value of ^ into the damping ratio comparator, and 
the adaptation servo continued to change damper gain at a 
rate proportional to the indicated ^ - error. 

It vjas found that during the adaptation process, the const¬ 
antly changing damper gain caused errors in computed 
and ^ due to the fact that system parameters v/ere altered 
appreciably during one sampling interval. And, the comput¬ 
ed values of ^ were caused to vary in a random manner and 
were generally 0,C5 to 0,1 below the true value of damping 


ratio. 

3) As pitch damper gain increased, the effect of the time rate 

of change of system parameters on the calculations 

was magnified due to the increased sy^'tom loop gain. 

4) ^’hen system motions damped to zero, the computer V7as 

not supplied ^^/ith sufficient Inforration to perform calcu- 
lations, 


It can be seen how these factors comb?Lned to causa erratic operation of the 
adaptation rncchanisiu. As an example, if a value of ^ =0^6 v;cre computed 
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ai;cl set ^ con.parator, no nn\; convergent values of ^ were 

genereitecl for t%’o seconds^ a steady ^ - error signal » 0.1 would be fed 
to the adaptation servo and it would increase damper gain at a constant 
rate during tlie tvjo second periodo This is il lustrated in Figure 29* v:here 
was increased linearly from time « 3,5 seconds to 5,5 secondsc At 5,5 
seconds a convergent value of ^ , exceeding 0.7, w^as obtained, set into the 
^ - comparator, and damper gain began a dov/mzard drift* Hov/ever, the drift 
of Kd was generally upward because the erroneous damping ratio solutions 
were usually smaller than 0,7, 

The problem of drifting damper gain v?as solved by an alteration of the 
computer output roiitine in which outpvit logic was re*“programmed to 
set ^ = 0*7 into the ^ comijarator if no convergent solution were obtained* 
This action effectively ’locked/’ the adaptation servo at its present set¬ 
ting, and damper gain was not allowed to vary until generation of a nev; 
convergent solution* The above process will be referred to as the ’^adap^ta- 
tion-servo Icclr* in the remainder of this report. An important benefit 
accruing from the use of the servo loc’: was the fact that the plant v;as 
allowed to settle to a linear system during each period of nonconvergent 
calculations and thus the accuracy of subsequent convergent calculations 
was improved. It is recalled that, when the plant is in motion^ the calcu¬ 
lations converge during one quarter period in each half system cycle, and 
from this, it would be expected that damper gain ^^ould be changed in quart¬ 
er period steps during the adaptation process. This is borne out in Figure 
30 i^hich is a time history of the adaptation process with the servo lock in 
operation. 

The run shown In Figure 30 is identical in all other respects, except 
the use of the servo lock, to the run illustrated In Figure 29 that Is, 
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f’lc'ht conMtLoii 1, m'’ a. t-s'- on-rvecvo :'ain I.,/. Tt is saen that 
pitch daiTiper gain increrscc- ir, stcj s, as predictcil, and Ic settled out 
close to the predicted v'-lae o’* K, = 0.3'v4 in 2.3 seconds. After three 
seconds nmning time, the airplane motion damped cut and no values of ^ 
were available until t = 5 sec. During this Intei'val damrei' gain remained 
constant at 7^ = At t - 5 sec., the airframe V7as perturbed by an iiu' 

pulse and the comj^uter begari generating values causing damper gain to be 
moved into the 7, - 0.34 to 0,35 region again. The time history of notions 

Cl 

resulting from the sercnd impulse also serves to illustrate the behavior of 
the adapted airframe which displays a danping ratio = 0,7. 

A series of runs wavS male with the servo lock in operation using eight 
different values of adaptation servo gain ranging front C.l to 5,0, It 
was found from this series that values of exceeding I.O catised an ex- 

o 

cessive overshoot in damper gain during the adaptatlcn process and that 
values of K less than 1.0 resulted in a smooth adaptation transition, but 
time required to fully adapt was approximately inversely proportional to 
servo gain. On the basis of these results, it was decided that « l.C 
was optimum, and a second series of runs was conducted to investigate the 
operation of the adaptation loop, with this gain setting, at each of the 
11 flight conditions. Table X is a sample listing of the computer output for 
the three degrees of freedom run performed for flight condition 7, from which 
the reader can see the stepwise manner in which adaptation is accomplished* 
Items of interest in this program v/ere the time required to adapt the air¬ 
frame from its basic un-augm.ented condition to damping ratio ~ 0,7, and the 
steady state v«alue of damper gain upon completion of adaptation. These values 
are listed in Table XT for each flight condition. Other items listed are the 
theoretically pre licted gain for adaptation and the basic airframe short per¬ 
iod damping ratio, Fgiurc 30 shows that the steady state damper gain exhibit¬ 
ed small amplitude fluctuations about a mean value during operation^ and it is 
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Figure 30 . Time History of Airplane & Control Surface Motion 
and Pitch Damper Gain During Adaptation Process V/ith Adapt¬ 
ation Servo Lock 

Initial K^= 0, Adaptation Servo Gain = 1. 
Flight Condition 7 - M = 1.8 © 25,000 ft. 
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Table X 
of Freedom 
Servo Lock 


Llstin 5 of Computer Output from a Three Degree 
Flight Simulation V^ith Adaptation Loop Using 
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Table XI - Comparison ol Stedcly State Damper Gain Generated by the Adapta¬ 
tion Syster-;, with Theoretical Gain Required for ^ 0.7 at 
Each Flight Condition, tfr-f fo ilirframc I}^mplng 

Ratio also Tabulare«l. 


Adaptation Servo Gain 1.0 


Flight 

Condition 

Basic 

{ 

Time to 
.\dapt 

Theoretical 

ICci 

Experimental 

Kd 

1 

I'l = 0,9 
10.000 ft 

0,440 

1.4 Sec. 

0. 156 

0.16 

2 

M =1.2 
10.000 ft 

0.177 

2.2 

0.303 

0.29 

y 

M = 1. 
10.000 ft 

C.15& 

2.1 

0.284 

0.28 

4 

M = 0.9 
25.GCC ft 

C. 244 

1.7 

0.2 09 

0.26 

5 

M = 1,2 
25,000 ft: 

C. 100 

2,5 

0.372 

0.36 

5 

M = 1.4 
25,000 ft 

0.105 

2.5 

0.384 

0.38 

7 

M = 1.8 
25.000 ft 

C.09G 

2.5 

0.344 

0.34 

8 

M = 0.9 
45,000 ft 

0.141 

2 .2 

0.405 

0.44 

9 

M = 1.2 
45.000 ft 

0.066 

3,1 

0.538 

0.56 

10 

M = 1.4 
45.000 ft 

0.067 

2.9 

0.503 

0.51 

11 

M = l.S 
45.000 ft 

0.058 

2.8 

0.467 

0.47 
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It is sacn that tlie taticn times cor.arc i^vuraily tvl^h ^bcst cl the 

MIT system, and as would be expected, these Clight conditions v.lth "'ower 
values of basic damping ratio exhibited longer adaptation times* The ex¬ 
perimental values of steady state damper gain couM be determined at best 
with two place accuracy, but even so, reasonable agreement with the theo¬ 
retical values was obtained* 

Considerations in the Practical Mechanization of an 
Airborne Digital Damping Ratio Computer 

Perhaps the best method of reviex^ing the finally ’^optimized” Adapt¬ 
ation process is to follow through the flov; diagram of Figure 31* which 
illustrates the prograi^i sequence to be followed by the proposed airborne 
damping ratio computer. The pi-irnary difference between the airborne com¬ 
puter and the ^solut ion programmed into the three degree of freedom 
flight simulation is the requirement for an analog to digital converter for 
pitch rate input signals and a digital to analog converter for output signals 

to the adaptation servo* It was found from the simulation program that sirnj le 

# • 0 

trapezoidal differentiation of 0 and 0 yielded sufficiently accurate values 
of 0 and 0 for use in the calculations* 

From a short preliminary survey of computer requirements it was deter¬ 
mined that a special purpose digital computer having a 130 x^ord storage 
capacity and a repertoire of 12 arithmetic and manipulating instructions 
would suffice for this application. And, if it is assumed that each in¬ 
struction requires 5 microseconds for execution (a representative state- 
of-the-art-estimation), the sequence from raw data input to generation of 
a new adaptation-servo command will require 12,600 microseconds, or in 
other words, the total cycle will have a maximum possible repetition rate 
of 83 times per second* Hox;ever, the cycle utilized In this study \vas re¬ 
peated only 40 tines/second, which leaves a comfortable execution'speed 
margin* 
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Fi 5 ure 31. Flow Diagram of the Proposed Damping 
Ratio Computer for Use in the Adaptation Loop 
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VIIi - 


It has been demonstrated that a simple, special-purpose digital com¬ 
puter can be utiMzed to determine the natural frequency and damping ratio 
of the short-period mode of the free longitudinal oscillations exhibited 
by a representative high performance fighter aircraft. And^ it has also 
been shown that the calculated values of damping ratio can be used in a 
closed loop manner to vary the gain of a self adaptive pitch damper which 
is capable of maintaining constant airframe short-period damping ratio over 
the full speed and altitude range of the vehicle. The system investigated 
is capable of performing the adaptation process^ from basic unaugraented 
airplane to ^ ^ 0.7, in a i-eriod of time ranging from 1.4 to 3.1 seconds 
depending on flight condition; and this speed of operation was shown to be 
similar to that of the most recent MIT flight-vehicle self-adaptive system. 

In conducting an evaluation program with the proposed mechanization, 
in conjunction with a longitudinal three degree of freedom flight simula¬ 
tion, it was found that measurements of pitch rate rather than pitch angle 
were necessary to extract the short period characteristics from the combin¬ 
ed phugoid and short period motions. First and second derivatives of pitch 

rate are then determined in the adaptation computer and sampled values of 

• • • •. • 

0, 0, and 0 set into the iterative routine for determination of short period 
damping ratio and natural frequency. The system requires no perturbation 
signals or limit cycling for operation of the identification procesSj how¬ 
ever, the system will adapt only while the airframe exhibits pitching motion. 
It should be emphasized that the mechanization investigated in this study 
is not a complete longitudinal control system, in that there is no pro¬ 
vision for external pilot-input signals. However, it was proposed in 
Chapter V that Corbin''s fon/ard-loop-gain computing system could be con¬ 
verted from analog to digital form, and used in conjunction with the present 
computer routine to enable the Introduction cf exteinal ontrol ipp»<tSo 
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O'-' i!!' . i‘ i ■" ■ r ■' 'rLfG 

c<M . \ ii t?:' . cfli-. I » 1 *-^ ir o study cf '‘e ol Jorciag 

the*, a r;VKyc^^^x}ce of ^ O seilos cf Iterat ive calc^ilaci • is It is felt 

th^t^ If the ca^cu^£:tlons cctld be n'.ade 100 peixent coriveigertthe 
speed ai'id accuracy cf the ada^tatioii process could be greatly enhanced. 
Other lc«artovs \ 7 hlch vould require investigation before a practical mecheni- 
nation could be accorplished are: 1, the effect of noise (atmospheric tur¬ 
bulence) Introduced into the forward loop, 2, the effects of control sur¬ 
face rate and ari^plitude limiting^, and 3, the influence of transducer in¬ 
accuracies and noise on the oj;eration of the adaptation loop. 
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Table of Air\>lc.Lie Phyc-’^cal 


cha vac l ex 1st 1c ? 


Airplarje ipcIss 1^110 slugs 




19j'/i-2 r.lug-ft. 

b 


129, 5.33 slug ft,“ 



142, 398 slug-ftd 

^xz 


8,219 3lug-ftd 

Wing 

area 

394,84 ftd' 

Wing 

span 

36,58 ft. 

Wing 

inean aerodynamic chord 

12,24 ft. 

Wing 

aspect ratio 

3.39 

T , 

aC 


13,95 slug-ft." 

'^e 


1,259 rad/sec. (Cruising; 
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DERIVATIVE 


A 

DERIVATIVE 

. Cl- fu TCI. .•J.AEi...tT oiVA. .AA ' 

r.CAi'EivLS 

^16 a 

ivoning n.Ou.ent/ ail^rcn deflection lltile irr.^crt- -^e 

in dynoiiilc r^tabilitye 

r 

Rolling moment/ rudder deflection cross ccupllng terr¬ 

or little inrj.ortanca In dynamir stability. 

I 

Rclllng morcnt/ rol- rate - viscous daim Irg In rc/1% 
balances ailex'on raoment to province n convStant lolli'ng 
velocity for a given aileron defl^ctioHj has si.iil} 
danc ing effect on lateral-directional spiral modscro 

S 

r 

Rolling TiiOc> 2 at^ yanking ve'^ocity - contributes slightly 
to Dutch iol\ vnoGt lii'j ox'tant cause of spiral divergence 
mode. 


Rcllirng moiianc/ sideslip angle - dehedral effect^ large 
negative value v/ill produce good spiral characteristics 
and poor D^itch roll characteristics« 

^10 

Rolling moment/ sideslip angular rate^ usually negligibba 
in magnitude. 

C 

in. 

"t 

Pitching iioment/ elevator deflection^ control derivative, 
little importance in dynamic stabilityo 

C 

ni 

q 

ritching moment/ pitch rate, viscous damping in pitch, 
large effect cn damping of longitudinal short period 
mode. 

c 

Pitching moment/ x'ate of change ct angle of attack^ due 
to aerodynamic time lag and aeroelastic effects^ contributes 
to damping of longitudinal short period modec 

c 

mQf 

Pitching mt^ricnt/ angle of attack - produces lestoiing moment, 
cstablishvss frequency of Icngltudinal short period rode. 

C 

m 

pitching K'Oi.ient' sideslip angle - aerodynamic cross coupling 
dei'lvative, usuany^ negligible in inagnitudec 

C 

«<5 r 

Yawing moi'tjnt/ rudder deflection - control dcriv-^Clvc^ little 
importarsce. in dynau ic stability. 

C 

n 

r 

Yawing moment/ yaw rate - viscous dnr.qrLng in yaw^ major contri 
butor to dimving in Dutch roll (lateral-directIona 1 shor»: perl 
cd mode^ * 
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n 

r 




Ve t - :.ro'‘’ynauUL cron:?; ccup-ln^ 

CutV f.b!-^‘^s slightly to dA:*ipin|> Ditch roTK 


C 




»>* anc/ ar’.jle - produce?? diraccloiial re» 

storing \x.\ 2 Tiiy establi<?h^r froq-jerr-y cj"' Dutch roll nodti* 


C . 


Ya<77.n;x no ent/ '?ir>esltij augu'^ar rate - contributes sHghtly 
to 'h\..plng of Dutch roHo 


C Yawirp; r.o“;ent/ aileron deflection - control cross coupling 

^ derivative^ little Inportnnce in dynamic stabilltVc 


C_ ‘"ide force/ sideslip ait^.gle - contributes to damping of Dutch 

>:ol^ 


Side force/ roll rate - acrodynajnic cross coupling dei'iva- 
~p tive^ ..icgrftule usually negligible. 

Side forCeV yaw rate « magnitude usually negligible, 

t 

r 


. Side force/ ^IdesIIi: nygu'’ar rate - magnitude ustially 

* ^ neglIglL le. 


/*! 

L» 


cr 


c 


^(X 


c 



Lift/ angle of att >c]' - effects frequency of 1 ongl^tudinal 
phugoid Mode, ha^ iviijor effects on aircraft vianeuverability• 

Lift/ rale of change of angle of attach - due to aerodynamic 
ti-ie lag and aeroelastic effects^ has negligible effect on 
dynamics, 


Lift/ ehjvator deflection ' iragrltudc negligible 
tlv.;nal aircraft, has effect on tailless aircraft 
vjith short tail length. 


for convan- 
and those 


C 


D 


AirpLane drag ^'•oeffIciant - primary cortilb^itji to damping of 
longitudinal ;hvgold osciM:\tlon, 



Dreg/ angle of attach ^ negligible effect onlongitudinal 
phugoid I'lode because this mode takes place ac constant ang'^.a 
of attack. 










p . 3 a a IV 


tQ ATljr. Cr 

The eq^jitlons r th^ lv-o Tat 11* i -- p-ess- 

tone; Force ~ I ' x ixceleration arsd 

;-loirei;t “ I'iOh.er.t '£ Tne^t ui x Angular Acce'^erac lo.j 
From classical rachaiJes, the equations of motion til a rigid I'vO'^y In free 
space no aeLodyn^wic or gravity forces actlrg re. 

= m(u -rv) 

/^y - m ( V i- ru -pw) 

Fi - i' 

L,=> pl)(^ ^ 

M- ^ly -f- (p<^ ~ (f> Ay 

N= rli i- - (^-trp) lyi P<^ (ly-Iy) lyy 

Note that syi'bols £.re .iei^ijei in figure U It is seeii rh.-c rhe ^ast 
t\/o terms in each of the first three equations are th: pioii-ct an angu¬ 
lar and a lii'car velociLy and thus represent centtipeta.^cceleiaticu'.s dir¬ 
ected along che respective ^xes. iha 1.-st fou^ ceris f -: l the last 

three equations represent nc^ents created by inertial cross couplings These 
tervuS arise frcri o. of soar::es Includiug gyrese op-’" ■. j rC' esslor.o torques 
















































erclal cour 1 ’i'- • t ^ 

pointed cut thct tliC '' cei^ .t) ^.oi 3*:f.-:0-S ^^eccu'i* ^ r . = .strous 

results, v.hen the a-‘'rMtu!03 . nl ‘^aier> of atrora.it r^otlcn v^c i (ref«. 

3), It Is recalTal tba aitcr-i£t |:orsussc?^ veiilc"'’ .'alih 

causes the [uocaccs oc riiertia an 1 to ->e cer.. .. This e I . c vaaties 

tvo Inertia *5 coup Hug teixus iron each equatxcu. i/ext^ the c oia* osckis oi 
force due to gravity, aud the qcro Avna' ilc forces ixcul ’.^r? eiits ara Iticl ideJ, 


P 5 -r- 

r)f^(j =: — / 77 ^ v 5 />? 0 ^ E^n^ins Thrust 


mv 


= pv^ ^ru cosO s'in0 


p j 

/y7vv = -yov/ - mp COS& (ZOs(p ^ ^iT 


pl)^ - C^r (ly-I^:) + rl )^2 ^ PC^^y~St> 


<^Jy-pr (l^ -/x) i-y-'^IxT- '~ ^Xo ^ 2_ 


Ki'P ^Xt-^Xe 


P z 

P<^ (I%~ ly)-+■ P IXi ~ 


’fhe I,, CO terms re; resei 
Xe e 

which when tiultli. lied bv a 
sccpic precessioral torque 
terms are significaiit In a 
that tba wingspat.’, b- is u*: 
the roll and yavj equations 


t ‘he anguhor nor ontu\r. of rotation, -ngitM i»nrcs^ 
yav or pitch angular '»'e]ociL>^ , •icdu’v.os a gyro- 
cheut the pitch yaw ax^'s rtT'prct-X'iily* These 
turbojet powered aircraft* It should be noted 
ed -ir 'w-r itu. t a)0' fo-. -. vrrh, .. i" ^ 0'''aats i;i 

vd)eieas '.ear. ee.rody?,\;-r ^ ^■' - used iz 





noi..e5:< 


f r 1 .. icv.' - r 
A :-j c ut. . , ■ !.- .ulo.i- ... • '- • 

CO^”^ OTiaaC o *'' X 1.3'.v ■* V.I 5 V, ;i‘a- 5 , A • l . - 

into nor^i c::et- 


U ^ V ^ cxicjoloLne PorwccrJ velocity 

\f V Sin S ^ VB 
V Sin (X V (X 
V ^ B 

Vex 


ex nd B ? aerooYiiau'lc iu*" 

incllnatloti q\ forvar! velocity vector to the. boov Eho 

complete equations of motion with the aerodycoeffi'“tents i-.n cxpattded 
form are showii in table I, The reader will see that the aerodynanic rc- 
tary dcrlva,tlv6 terms are multiplied by cci^ain add\ttori.^il flight Ooodltt.wia 
parnrneters in order tha‘<r the v?:liie of the derl\.-^t Ivi rem;atn as ce^'' 

as possible over a wl'le range of flight conditions* exafnpls l9 ^ he use 

of the term ^ with derivatives taken with respect to roll "Ifva.lar 

terms for other deriv..itivas ^ra listed in thv^ table of sv^^iboSs^ 

One final sec of ktiquations most be used to conv-^^rsiei: of 

angular rates about the body axes to ?mgolar rates wifh t'eSi'-iet tO' the 
earth coordinate r>ysrc,mo The necessary 3uiet equ.^.t tons are 5.)- 
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0= CO ^ cf) — r sin (/) 

0 - p T Ks}n<^ tan <p ^ r co^ (p tan d 
^ sJn ^ r cos 0) S(S(C 9 

• » 

G and 0 can be integrated to obtain 0 and 0 v^hich ai^e used in the equa¬ 
tions of ntotlon. 
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DE^IVATIOi] OF T'E LO.%GI'FiJ!;I.-AL 3HORT-V£IlToi.* 


^.H S FE F'^.‘. iC 110., 


inr 


EOUiiTIONS OF HOIION 

Ou the two lon^ltacUiir: 1 modes of airplane mot ton, the phjgoid and 
the short-period, only the latter mode Is of interest in stahillcy augr«jeritn» 
tion considerations because the pbugold nornjslly does not cause any deter-^ 
ioratlon of handling qualities^ In Section III, a transfer function of the 
mode of interest in stability augmentation syst^3ms must be included in the 
bloc:k diagram of the proposed [tich damper to enable an analysis of the 
systenu Thus, the longitudinal sbort^periocl transfer function must be de¬ 
rived from the airplane equations of motion. The Jegrees of freedon* of 
interest being pitcli and noruial translatioti for which tlte corresponllng 
equations are shox/n helov;. 


cc^e ^ cosQ - cx 


The pitch da'cping derivatives^ snd us.r-^lly i-jdda-i together 

to fomi ai. equivalent engine torque and tail, lift cern-s 

have been neglected because of their Tot' ii'rigic’jdes; rissi.iv” {.'’c-ns rieces- 
sary for the derivation of the transfer rur.ction are: T, :-r- oqui 1 tbriu-.n 
condition exists a-ong the ' axis vhen airplane weight ^ airf lift, 

W=- 1. 

The difference between CC and CC -trim vLll be denoted cr, 
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r 


Ofia.U. ions of n‘o{'Lon then hecore: 


0 = 


he 


Crn -. ‘t 


A 

ZV/y 




-h 


y)c 

Ly 




Oe = 0 - 


A 


nrV 


C,CL. 


CC 


he next step is the IntrovluctIon of the liffcrentlal operator^ S, and the 

0 


elimination of angle of attack from the equnti^ons in order that the 
equation may be obtained. 




c ^ 1. / Ac^ 


96 CC, 


C03 - QS lYj^ 


or: 


CC, == 


9S 


S -h 


^L, 


nov 


Substituting CC, 


i Ac^yr. Ac €6 




Ac ^ : 

Zy ^ 


e 


The — transfer function Is: 

Ao 


0 


i-y Cm-,^{Sy- ■ynT'^La.) 


r / oZ ^ r ^ p - Ac ^ AAit-r A ^ Ao \ 

ZViy ^fO^O ^Lec^ lAT^^^cc) 




Rearranging 


0 


Ac 

/v 


Cmi^ (^S 


A ^ 

^LocJ 


mv 




TIB 




























h c> 


*s ‘-^hiv n hc-'V 


' 1 ' 


0 




*■ 'n • ' ' ' t: r,_<■’• 

/^e(SrTe) _ 

6(6- t2^uJn6 -rOUr,^) 

*h jTxS: 


•'"he: n Li-i ■’•r.ne Ci'*evnr,ur e Tfect LVt^nesj-’ 
7e the night path tu.ie constant 

^ i s the sh o r (: - r e r 1 o darnIng ra t ? o 


Is the short “per i*o4 natural frequency 


t can be seen that the steady 


state transfer function is 


assent Lnlly 


pitch rate coimKind systen^ 


for a given horii^ontal si 


zev -leflectIon, 
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.1C EVSTEM 


:i‘rnA- 


A'’thou.;h nnoJo;^ cor^. uter? ^.^eneraVy used to sln.ulate syst3ui 

dynamics, the systems i.nvest l^nteo In this study v-ere s.it^iulated oit a stored 
progratr?, gonetal purpose digital computer using techniques similar to ana¬ 
log computer programs. The techniques arc similar in that both methods use 
sone form of Integration coiub'nod vjith other operations such as summation, 
subtract ion, and multi]pl lent Lo^ to sedve a set of slTiultaneous differential 
equations of motion; each equation of motion representing a degree of free¬ 
dom of tlte system* An analog Cvmputer program may he regarded as a contlnu- 
oxis laraHel solution of the equations of i-otlon, x>;hereas the digital sjethod 
employs a step^ by step serial solutloiu A number of digital integration 
techniques are available, ranging fron» simple rectangular integration to 
highly accurate, biit time consul ir.g r.ethods, such as the Runge-Kutta numeri¬ 
cal integi-ation. The tcchnlqu.e selected for this study is a modified form 
of the trapezoidal rule and Is shoivn hclov;, 

!<■.- X.-: *i-[ •3 

u-here: n-h denotes values computed in previous cycles 

h is the v/idth, in time, of the coi-iputation interval, 

A major difference betveen the digital solution and the concept of the 
analog flow dlagratn lies in the fact that the digital method perforins in¬ 
tegration of lower ordered derivatives first. As an example, the steps 
in the computation cycle of a seconl ordv^r system vSoUition are shown he'^ow: 


" X/?-/ 

X/7 - X/?-/ ^ LJO-/ ” X/)_2*7 


'20 







The cot 1 Ion r . cV. f.5? v'ej'eated at ^ high rate and tho values of X^XjS)( 
genoratod In each rcr^iresene lolnts on a time history of syston motion. 

^Jhen a more ccnivlex systeu: having a greater numhcr ol euuations of motions 
is simulated, the Integrations for each equation of motion arc perforrued 
in a serial manner such that one ]ass is executed through all equations 
before recycling. 

The accuracy of this technique is, of course, a function of the 
vidth of time incren>ent betv.een comjuting cycles, and it v;as noted that 
when using t1»e method described above, integrator errors were manifested 
as an undarnping of the system dynamics. In some lightly damped cases it 
was found that cumulative integrator errors actually caused the solution to 
diverge, the corrective action being the use of finer calculation inter¬ 
vals. It has also bean noted that the required fineness of calculation 
intervals is a function of both the period of the system mode of highest 
frequency and the damping of this mode. A rough rule of thumb is the use 
of 20G divisions per period for ^ = 0.1, dovTn to 50 divisions per period for 
^ = .5, to obtain one percent solution accuracy. 

The following five pages are sample listings of two of the FORTllAIl- 
programming-language programs v;rittcn for this study* The first is a 
second order systei^. and the second a three degree of freedom flight simula¬ 
tion. The programs V7ere run on both the GDC 1G04 and the IBM 7090 general 
purpose digital computers. 
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SECOND ORDER PLANT IDENTIFICATION ROUTINE 5 


M=i 

-1 .R£AD. II'iPUT TAPE 3 
1A,STHETA, ALPHAT, 

2 FORMAT (F4.i, 16, 

read input tape 3 

3 . FORMAT ( 6FC-.3, F8 

WRITE OUTPUT TAPE 

-4—EORMA T_t,l Tiiy / 5 8H 

1NUMBER= F4.1, i2H 


ALT, DENSTY, VELX, DYNPRS, LUEJA, C TUET 


-INITLAL- 

J = 1 

-Ji = l 

K = 1 


CONDITIONS 


,2, ACHM, 

THRUST 

F8.7, F7.2, F8.1 , 4F7.6, F8.2) 

CMIT, CMUA, CMA, CLA, CLDELE, CDU, TIIEDO 

• 4 ) 

2,4, ACHM, ALT 

-THREE. DEGREE OF-. FREEDOM- FLIGHT SIMULATION MACH 
ALTITUDE= 16) 


Ll=l 

-- TIME=0. 

SRVGN=.25 
. —DAMPR=.7. 
DKPGN=0. 

-DAMPR2=.4_ 

ACCLN=0. 

- OMEGAN=0. 

50 ALPHA=.5 

_ALPHD=0._ 

ALPHD1=0. 

_ THETD.= 0.__ 

TH£TD1=0. 

- THEDDi = THEUD 

THTTD=0. 

_THTTD2 = 0. . 

.VELCTY = VELX 

- ^VELXD=0^-J _ 

VELXDi=0. 

_ACTP0S=0,_ 

ACTVL=0. 

_ACTVLi = 0. 

ACTCL=C. 

_ AC.TCLi.= 0.. _ 

SIGNAL=0, 

ALPHA2 = .5_ 

ALPHD2=0. 

._THETA2 = . 5 

THETD2=0. 


_THEDD2 = THEDD _ __ 

5 WRITE OUTPUT TAPE 2,6 7'“ “ -- 

„6 FORMAT (105H TIME-SEC _ - THETA THE.TADOT FWD VELOCITY ACT 

lUATQR PCS NORM ACCEL DAMPER GAIN OMEGAN OAMPR) 

7 WRITE OUTPUT TAPE 2,8,TIME, THETA, THETO,VELCTY, ACTPOS, ACCLN, 
IDMPGN, OMEGAN, DAMPR 

8 FORMAT (9F12.''6) _ 

IF (L-iOOO ) 9, 12.» 12 ' ’ ' ' - - 

9 L=L+i 

IF{Ll-200151,52,52 

51 L1=L1+1 ' 

GO TO 14 
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SECOND ORDcR PLANT IDENT IFICA1I ON ROUTINE 3 


!d2 L1=1 

- . GQ__T-0-^5 0- _ . __ _ _ ___ 

12 CALL EXIT 

ACTUATOR SIMULATION 

14 ACrPas = ACTPOS+.002 3*(3.*ACT\/L-ACTVLi ) 

ACTVLi=ACTVL 

ACTVL = ACTVL+. 002 5<^ ( 3.-»ACTCL-ACTCLi ) 

_ACTCLi = AC.TCI_ . , . __ 

ACTCL = -ACTVL»49.-ACT POS*12Z5. + SI GNAL<-122 5. 

INTEGRATION OF LON ITUDINAL ACCELERATION TO OBTAIN \/ELjCITY 
VELX = VELX+.0025* ( 3 .-^VELXD-VELXDl ) 

_VELXDI=VELXD 

VELXD= (THRUST/illO. ) - ( AOVbKM»COO ) - ( AOVERM* ( ( CL A« ALPHA) <^*2. ) / ( 3,14. 

_I59j:3.-\9) )rJ3Z..i8»SINF( THETA) ) ____ __ 

VcLCTY=VELX 

. -.-CALCULATION OF AERODYNAMIC PARAMETERS DEPcNDENT ON VELOCITY 
DYN.PRS = DENSTY»( VELCTY*»2. ) «-394.34 

_ ACIYY=DYNPRS«12.24/129533. 

AC2VIY=(ACIYY»12.24)/(2.«VELCTY) 

_AMV=DYNPRS/(V£LCTY«1110.) __ 

AGVERM=DYNPRS/1110. 

SOLUTION OF EQUATIONS OF MOTION 
THCTA=THETA+.0025*(3.*THETD-THErDi) 

_THETD1 = T.H£TD _ _ __ . 

THETD=THETD+.002 5*(3. i^THEDD-THEDDi ) 

_THEDD1 = THEDD-_ ... _ _ _ _ „ ... 

THEDD=(ACIYY*CMIT*ACTPOS) + (AC2VIY»CMQA*THETD) + (AC IYY^CMA*(ALPHA-AL 
__-lPHAT ) ) 

ALPHA = ALPHA+.0025‘^ ( 5. *fALPHD-ALPHDl ) 

ALPHOl = ALPHD 

.. ALPHD=THETD+ ( 32.1 8*COSF ( TH £TA ) /VELCT Y ) - ( AMV*CLA<t AL PHA ) 

__ CALCULATE NORMAL ACCELERATION . __ _ - ___ 

ACCLN=(CLA*ALPHA»DYNPRS)/(ii10.*32.18) 

_ PITCH DAMPER SIMULATION , _ _ 

SIGNAL=THETO*OMPGN . . • • - 

_iTIME = TIME + .005 v.> 

DIFFERENTIATE PITCHING ACCELERATION 

_THTTD=(THEDD-THEODl)/.005 .. ... 

ADAPTATION SERVO 

DMPGM=DMPGN+( (.7-DAMPR)* SRVGN*.005) 

' IF(J-5)i5,i6,16 • • 

15 J=J+1 

GO TO 14 

16_J = 1 .. . ... . _ _ 

CALCULATE NATURAL FREQUENCY AND DAMPING RATIO- ASSUMED DAMPR=.4 
DAMPRi=DAMPR2 ‘v 

21 GAMMA=DAMPR1*THEDD2/THETD2 

BETA=SQRTF(GAMMA**2.-THTTD2/THETD2) 

OMEGAN=-GAMMA+b£TA 

DAMPR1=- ( THT'i IJ/( 2. >-QMlGAN*THEDD) ) - ( ( OMEGAN*THETO ) / (2.*THEDlj) ) 

27 IF( Jl-20)20,3‘1,3i 

30 J1=J1+I 
GO TO 2l 

31 Ji=l 

I F(OMCGAN-8. )25,26,25 
2 5 IF(nMCGAN)26,52,3_2 
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SECOND ORDER PLANT IDENTIFICATION ROUTINE 3 

I 

32 IF(DAMPRl-l.)33 t26,26 

_33_iFLDAMP.Rl ) 26,.o4,34 _ _ _ __ 

34 DAMPR = DAMPR.l "■ 

DAKiPR2 = DANPRl ^ ___ 

GO TO 29 

26 DAMPK1=DAMPR2 __ 

35 GAKMA=DAMPKi*THEDD2/THETD2 

.BETA=SQRTF{GAM(^A**2,-THTT02/THET02) _ _ 

OMEGAN=-GAHMA-c>ETA 

DAMPRl=-( TIiTTO/ (2. ^fOMEGAN*THEDD) ) - ( { OMEGAN*THETD ) / { 2. «THEDD ).)_ 

36 I,F( 01-20)37,38,38 

37 J1=J1 + 1 •“ - _ _ __ 

GO TO 35 

.38_,J1 = 1 _ ___ _ 

IF(OMEGAN-8.)39,40,40 

39 IF(0MEGAN)40,4l,41 . __ 

41 LF(DAMPRl-i.)42,40,40 

.42 IF(DAMPR1)40,54,34 _ _I_ 

40 DAMPR=.7 

,_WRITE_OUTPUT_ TAPE 2,43 _ __ 

43 FORMAT (31H NO SOLUTION FOR DAMPING RATIO) 

29 THETD2=THETD . . 

TM£DD2=THED0 

THTTD2 = THTTD .. 

IF(K-iO)17,18,18 

GO TO 7 

-18 i< = i ______ 

GO TO 5 

END(i,0,0,0,0,0,0,0,0,0,0,0,0,0,0). .... 
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PLANT IDENTIFICATION ROUT INE_SECOND ORDER 2 


\ 


XEC 
1 = 1 
J=1 
K = 1 

XD0=0. 

XCDO=G. 

XC = 1. 

XDN=0. 

XK=1,_ _ 

TIKE=0. 

XDOO=e« 

XDDOO=0. 

l^RITE OUTPUT TAPE 2,1 
FORMAT (1H1//41|H PLANT 
CAMPR=.4 

XDDN=-2.4»XON-36,*XN 


identification ROUTINE - SECOND ORDER) 


XDINCR=(XODN+XCDO)*.00025 
XDN=XCO+XDINCR 
XCDO=XCDN _ 

XINCR=(XDO+XDN)*.0C025 
. XN=XN+XINCR 
XDO=XDN 

TIME=TIME + .Ci005 „ 

IF(1-1000)4,5,5 

4 1=1+1 _ 

GO TO 2 

5 1=1 _ _ 

VvRITE OUTPUT TAPE 2,6 

-6 FORMAT .( 1H //55H 
IEMENT) 

-_ ■ 

ITERATIVE CALCULATION OF NAJ.UR.AL_FReQUENCY AND DAMPING RATIO riVF 

A STARTING VALUE OF DAMPING RATIO = 0.4-.RATIO_GiyE 

WRITE OUTPUT TAPE 2\8 _ ‘ . 

8 FORMAT (25H CMEGA-N 

9 ALPHA=CAMPR*XDOO/XO 


TIME_ 


ACC ELER A11OJM V E LOC I T V_01 S P L A 


DAMPRAT) 


*XDN)) 


BETA=SeRTF{ALPHA«»2.-XDDOO/XO) 

0MEGAN=-ALPHA+8ETA _ 

DAMPR=-(XDDN/(2.»0MFGAN*XDN))-((0MEGAN*XN)/(2 

OMEG=ABSF(OMEGAN) _ _ - 

IF(OM£G-100, ) 15,16,16 ^.. . 

15 DMP=ABSF(DAMPR) 

IF(DMP-1.)17,16,16 -- 

16 WRITE OUTPUT TAPE 2,18 

18 FORMAT (19H SOLUTION DIVERGING) ^ 

DAMPR = .4 _ _ _ _ • 

20 ALPHA=CAMPR«XDOO/XO 

BETA = SCRTF(A'LPHA«*2.-XDD00/X0) 

0MEGAN=-ALPHA-BETA_ _ 

Oi5eG = IsSR 

IF(0MEG-100.l21,22,'>2 

21 DMP=ABSF{DAMPR)_ _ 

IF(DMP-1.)23,22,22 “ ~--— - 
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PLANT IDENTIFICATION ROUTINE SECOND ORDER '2 

22 V^RITE OUTPUT TAPE 2t24_ 

24 FORNAT (19H ST)LUTICN DIVERGING) 

GO TO 12 __:_ 

23 WRITE OUTPUT TAPE 2,25, OMEGAN, DAMPR 

25 FORMAT (2F13,6) . . _ _ 

IF(J-20)26,12,12 

26 J = J+1 ___ 

GO TO 20 

17 WRITE OUTPUTTAP E _2, J 0., _0 MEG A Nj„D A M.P.R _ 

10 FORMAT (2F13;6) 

IF(J-2C)11,12,12_ 

11 J = J+1 ■ 

GO -TO 9 _ r: _ 

12 J'=l . 

XO=X'N _:_ 

XOOO=XCN 

XCDOO=XDDN _ _ 

IF/(K-40) 13,14,14 

13 K = K+1 __ 

GO TO 2 

14 CALL EXIT __ 

END( 1 ,0,0,0,0,0,C,0,0,0,0»u,0,i«>,v) 
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